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AB STRACT 

This Is the Third Semi-Annual Status Report on Grant NAG 3-74 and 
discussion Is on those activities undertaken during the first naif of the 
second year of the study. I^asurements and computations are being applied 
to an axlsymmetric swirling flow, merging from swirl vanes at angle $, 
entering a large chamber test section via a sudden expansion of various 
side-wall angles a. New features of the present year's study are: the 

turbulence measurements are being performed on swirling as well as non- 
swirling flow; and all measurements and computations are also being performed 
o a confined jet flowfield with realistic downstream blockage. Recent 
activity in the research program falls into three categories: 

1. Time-mean flowfield characterization by five-hole pitot probe 
measurements and by flow visial ization. 

2. Turbulence measurements by a variety of single- and multi -wire 
hot-wire probe techniques. 

3. Flowfield computations using the computer code developed during 
the previous year's research program. 

Section 1 of this document describes this recent activity in the above 
subsections of the study, while Section 2 gives a brief summary of the 
present status and a list of recent publications. Appendices A through E 
contain a recent report, an M.S. Thesis, and three research papers. 

Finally Appendix F contains a report on the first six months of the sub- 
contract activity av. Dynamics Technology. 
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1. RECENT PROGRESS 


1 . 1 Flow Visualization an d Time-Mean Ve lociti es 

Time-mean flowfield measurements have been performed with a traversing 
five-hole pitot probe in a confined jet with downstream blockage. These 
measurements have been performed on a nonswirling confined jet and are cur- 
rently being done on various swirling flows. It was determined that the 
downstream hiockage had a significant effect on comer recirculation patterns 
as the blockage position is closer and closer to the inlet. The data obtained 
from these experiments will be Important in evaluating the predictive capa- 
bility of the k-c turbulence model and STARPIC computer code flowfield 
calculations. They will also assist in turbulence model developments. 

As part of this activity, careful measurements of the flowfield immediately 
downstream of the swirl pack have been made. These measurements have been 
performed and documented for various blade settings. Briefly, the measurements 
show that the flowfield typically has a swirl flow angle which is a ftw degrees 
smaller than the blade angle. Also, the measurements show a larger axial flow 
in the outer annular portion of the duct; presumably "thrown" out by centrifugal 
forces associated with the swirl as the flow passes through the swirler. 

Reports on tnese activities will be amalgamated later. 

Regarding flow visualization, photography of neutral ly-buoyant helium- 
filled soap bubbles, tufts, and Injected smoke helps to characterize the time- 
mean streamlines, recirculation zones and regions of highly turbulent flow. 

1 . 2 Hot-Wire Anemometry 

A significant effort has been made in the application of two hot-wire 
methods to the measurement of time-mean and turbulence properties, and a third 
method is being readied. The experiments being performed are designed to 
provide the information necessary for turbulence modeling development in the 
confined jet facility. 
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1. On e-W ir e Metho d. Single normal hot-wire measurements in non- 


swirling flow have established that the experimental technique is pro- 
ducing reliable results. Six-orientation hot-wire measurements have 
been made from which a method has been developed to make estimates of 
all components of the time-mean velocity vector and the Reynolds stress 
tensor. A major portion of these measurements are now complete for the 
nonswirling flowfield, and corresponding measurements in swirling flow- 
fields are currently in progress. The effect, of downstream blockage 
is also to be investigated under nonswirling and swirling conditions. 

Currently, the nonswirling case without downstream blockage has been 
extensively evaluated, and the results of this activity are reported in 
an MS Thesis, which is included as Appendix A of this report. 

2. Two-Wire Method . Although the single hot-wire (six-orientation) 
measurements are convenient and provide a great deal of information on the 
time-mean velocity and on the kinetic energy of turbulence, they have a 
shortcoming in that measurements of the Reynolds shear stresses are less 
accurate than desirable. Conseguently , crossed hot-wire capability has 
been developed to measure turbulence properties. This technique has been 
used for the nonswirling flow without downstream blockage, and is soon to 
be used with downstream blockage. A new cross-wire probe is on order which 
will enable a more accurate evaluation of the turbulent shear stress than 
currently possible. The results of our present activity is included in the 
report in Appendix B. 

At the moment it seems that application of the crossed-wire method to 
swirling flowfields will not produce measurements of all the turbulent stresses 
unless the probe is rotated through multiple orientations. However, the 
six-orientation single-wire method will yield appropriate estimates. Presently 



this technique is being used to measure turbjlent stresses in a confined 
swirling flowfield, and the results will be reported in the next Semi-Annual 
Status Report. 

3. Three-Wire Method. One of the major results of the six-orientation 
single hot-wire technique was to show that estimates of the turbulent shear 
stresses typically have large uncertainties. In fact, the uncertainties 
are probably so large as to preclude a definitive evaluation of turbulence 
models used in the computations. Believing this to be the case, a series of 
hot-wire measurements have been planned, using a three-wire, hot-wire probe 
with direct computer interface and data reduction. These measurements are 
to be performed by subcontractor Dynamics Technology, Inc. under the direction 
of Dr. Dennis K. McLaughlin. Me is being assisted by Salim I. Janjua who 
joined the company in January 1982 after completing his MS degree at Oklahoma 
State University. Dynamics Technology has been assembling a swirling jet 
facility which will duplicate the flowfield produced by the Oklahoma State 
University facility. However, it will be instrumented differently, with the 
three-wire hot-wire probe and its three-degree of freedom probe drive being 
major constituents. Adequate progress has been made on the facility assembly 
in order to complete the experiments during the first half of 1982. A 
report on these activities appears in Appendix F. 



1 . 3 Computer P rediction 

An advanced computer code Is used to predict corresponding confined jet 
flows. The calculation method includes a stairstep boundary representation 
of the expansion flow and a conventional k-e turbulence model. The predictions 
include recirculation zone characterization and mean streamline patterns, 
which are being compared with on-going experimental studies. The general 
development was discussed in the previous Semi-Annual Status Report, and the 
details are included In the Ph.D. Thesis of D. L. Rhode and its associated 
report, NASA CR-3442. Currently, predictions are being made with the standard 
k-c turbulence model and no downstream blockage. These constraints will be 
removed later. Predictions now include realistic inlet conditions. 

An indication of predictive capability is shown In Fig. 1, taking 
measured u, v and w profile* very close to the 10-blado variable angle swirier 
located 4 cm upstream of the expansion x/D=0 location. They may be compared 
with measurements shown in Fig. 2 [although the experimental data shown at 
x/D=0 are actually 4 cm downstream of the swirier, and cannot therefore be 
compared directly with the inlrt profiles used in the code and shown in Fig. 1], 
The corner recirculation zone is predicted to be slightly too large, and the 
central recirculation zone slightly too narrow. However, it must be remarked 
that worse predictions occur it the inlet radial velocity v is assumed to be 
zero, or flat inlet profiles, or solid body rotation inflow, is assumed. The 
present predictions are superior to previous ones by Rhode, who used flat inlet 
profiles initially, then realistic u and w profiles but with v equal to zero. 
There is a need for precise inlet flow data to be recorded in experimental work, 
and clearly turbulence model developments are needed for realistic flow simu- 
lation. Similar computations have been generated for other swirl cases, and 
a complete discussion on these activities is to be prepared. 
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Oblique to the present research, a simple technique for the solution 
of transient 2-D and 3-D flowflelds has been computerized. The k-e turbulence 
model has recently been added, and species diffusion and buoyancy are included. 
A paper on the application of the 2-D code to vortex controlled diffusers Is 
given In Appendix C, while Appendix D Illustrates two applications of the 3-D 
turbulent code. Configurations Investigated are local destratification near 
the release structure of a reservoir and deflection of a jet entering normally 
into a uniform cress-flow, showing good agreement with experimental data In 
both cases. 
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2. PRESENT STATUS 


2 . 1 Sum mary 

There are three major components In the present research program: 

1) time-mean flow characterization, 2) turbulence measurements, and 3) flow- 
field computations. Numerous activities are underway in all three components 
and results are already being produced and docunented. No major difficulties 
are anticipated in meeting the goals and objectives of the NASA grant and 
making a substantial contribution to the understanding and prediction capa- 
bility of complex turbulent swirling flows. 

2 . 2 Publi cations 

The following list of publications covers research conducted with the 
support of the grant, and does not include those listed in previous Semi- 
Annual Reports: 

1. Lilley, D. G., Computer Modeling in Ramjet Combustors. AIAA 
Journal, Vol . 19, No. 12, Dec. 1981, pp. 1562-1563. 

2. Janjua, S. I., Turbulence Measurements in a Complex Flow- 
field Using a Six-Orientation Hot-W’re Probe Technique. 

M.S. Thesis, Oklahoma State University, Stillwate . OK, 

Dec. 1981. 

3. Busnaina, A. A. and Lilley, D. G., A Simple Mnite Difference 
Procedure for the Vortex Controlled Diffuser. Paper No. AIAA- 
82-0109, Orlando, Florida, Jan. 11-14, 1982. 

4. Busnaina, A. A. and Lilley, D. G., A Basic Code for the 
Prediction of Transient Three-Dimensional Turbulent Flow- 
fields. Paper No. ASME 82-FE27, for presentation at the 
AIAA/ASME 3rd Joint Thermojhysics, Fluids, Plasma and Heat 
Transfer Conf. to be held in St. Louis, Missouri, June 7-11, 1982. 
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FIGURE LIST 


Fig. 1 Prediction of axial u and swirl w velocities for the sudden 
expansion «*90 deg flowfield with swirl vane angle 4>=45 deg, 
using measured 1n*®t values of u, v and w. 

Fig. 2 Five-hole pitot probe measurements for the corresponding case 
of Fig. 1. 
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CHAPTER I 
INTRODUCTION 

1.1 The Combustor Flowfleld Investigation 

Understanding the fluid dynamics of the flow In a gas turbine com- 
bustion chamber has been of great concern to designers In recent years. 

A gas turbine combustor, shown Figure 1, Appendix B, must burn fuel 
completely, cause little pressure drop, produce gases of nearly uniform 
temperature, occupy small volume, and maintain stable combustion over a 
wide range of operating conditions. The designer has a formidable prob- 
lem In aerothermochemistry, and more thorough and accurate procedures 
can help In accomplishing the design objectives more quickly and less 
expensively in the near future. 

Intensive research Is being carried out at Oklahoma State Univer- 
sity on the si.uject of gas turbine flowfleld Investigations ir the 
absence of combustion. Figure 2, Appendix B, shows the characteristics 
of the simplified flowfield being Investigated. Flow outers through a 
jet c* diameter d into a tube of diameter D, after being expanded 
through ar. a .-*q 1 e a. Before entering the tube, the flow may be swirled 
by a swirler located upstream of the Inlet plane. The flowfleld Is 
presently being investigated using various methods of approach, such as 
computer modeling of the flowfleld and flow visualization for both swirl 
and nonswirl conditions (2, 3). 
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1.2 Previous Experimental Studies 
on Expansion Flows 

Several studies on time-mean flowflelds of the type just described 
have been carried out using various turbulence measuring techniques (4- 
11). Unfortunately, most of the techniques used do not give complete and 
detailed Information about the flow In terms of all Its time-mean and 
turbulence quantities. There Is a strong need to obtain all the turbu- 
lence quantities In a complex flowfleld using a minimum amount of instru- 
mentation and without causing a great deal of Interference with the flow. 

1.3 The Turbulence Measurement Problem 

Turbulence measurement in a complex flowfleld has always been a 
complicated problem encountered by engineers. In the past, turbulence 
phenomena have been discussed by various authots in detail and various 
methods of turbulence measurement have been suggested (12-15). One of 
the most widely used Instruments to obtain turbulence quantities is the 
not-wire anemometer. The most common of all hot-wire anemometers is a 
single hot-wire. When used at a single orientation and in a two- 
dimensional flow, a single hot-wire can measure the streamwise components 
of the time-mean velocity and the root-mean-square velocity fluctuation 
at a particular location in the flowfield. A two-wire probe can be used 
to determine the time-mean velocities, streamwise and cross stream tur- 
bulence Intensities, and the cross correlation between the two components 
of the velocity fluctuations (16-18). To measure the three velocities 
and their corresponding fluctuating components in a three-dimensional 
flowfields such as encountered in combustor simulators, there are two 
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methods that can be employed at a point In the flowfleld: 

1. A multi-wire probe used with a single orientation. 

2. A single-wire probe used with a multi-orientation. 

Multi-wire techniques, with three hot-wires mounted on the same base 

so that they all lie within the same volume of the flowfield, permit the 
necessary three sets of readings to be made simultaneously. The require- 
ment Is to determine all three components simultaneously. The main dis- 
advantages of such a technique are: 

i. It requires three closely matched anemometer units. 

11. The probes Inter? »re with each other unless they are 

carefully placed relative to the time-mean velocity vector. 

ill. The spatial resolution Is poor because of the large size 
of the probe assembly. 

1v. Heat can be coi.vected from one wire to another giving 
biased readings. 

Multi-orientation of a single hot-wire Is a novel way to measure the 
three components of a velocity vector and their fluctuating components. 

A method devised by Dvorak and Gyred (19) uses a single normal hot-wire 
oriented at three different positions such that the center one is 
separated by 45 degrees from the other two. The velocity vector at a 
location Is related to the three orthogonal components using pitch and 
yaw factors as defined by Jorgensen (20). The data are obtained in the 
form of mean and root-mean-square voltages at each orientation. However, 
the measurements done with a single wire do not supply all the infor- 
mation needed to obtain the turbulence quantities. Therefore in addition 
to a single wire, Dvorak and Syred used a cross-wire probe to obtain the 
covariances between the voltages obtained at adjacent hot-wire 
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orientations. A cross-wire probe, two wires mounted on the same base 
and separated by 45 degrees from each other, poses the same problems as 
already discussed for a multi-wire probe. 

King (21) modified the technique developed by Dvorak and Syred. 

His method calls for a normal hot-wire to be oriented through six 
different positions, each orientation separated by 30 degrees from the 
adjacent one. Thus, one measures mean and root-mean-square voltages at 
each orientation. The data reduction is done usi.*g some assumptions 
regarding the statistical nature of turbulence, making it possible to 
solve for the three t1m*;-mean velocities, the three normal turbulent 
stresses, and the three turbulent sheer stresses. Having obtained these 
quantities, one can in addition calculate the kinetic energy of turbu- 
lence. Various recent studies discuss the turbulence measurement prob- 
lem, with emphasis on hot-wire and laser anemometer applications to 
swirl flows (22-23). 


1.4 The Scope of the Present Study 

In the present study, the six-orientation single normal hot-wire 
technique Is being employed to obtain the turbulence quantities in the 
combustor simulator confined jet flowfield. Measurements have been 
carried out for nonswirling flow with expansion angles of 90 degrees 
(sudden expansion) and 45 degrees (gradual expansion). 

Chapter II gives background Information on the various components of 
the exper mental facility and the instruments employed for the hot-wire 
measurements. 

The response equations using King's approach are given in Chapter 
III. Certain deviations from the procedures suggested by King are also 
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included in this chapter. A thorough uncertainty analysis of the techni- 
que is carried out in order to judge the accuracy and the reliability of 
the slx-oi ientation hot-wire technique. The salient features of the 
analysis are discussed in Chapter IV. 

Turbulence quantities obtained, using this hot-wire technique, are 
part of Chapter V which discusses the results In detail. Some of the 
turbulence quantities are compared with measurements done by Chaturvedi 
(5) using cross-wire probe in a corresponding flow situation. Chapter VI 
concludes by summarizing the major achievements of the present study and 
suggesting some avenues for further research activity. 


CHAPTER II 


EXPERIMENTAL FACILITY AND INSTRUMENTATION 
2.1 Idealized Flowfield 

The facility, designed and built at Oklahoma State University, is 
a simulation of a typical axisymmetric combustion chamber of a gas 
turbine engine shown in Figure 1, Appendix B. The schemat.r of the test 
facility with idealized flowfield is shown in Figure 3, Appendix B. 
Ambient air enters the low-speed wind tunnel through a rubber foam air 
filter. Next the air flows through an axial flow fan driven by a S h.p. 
varidrive motor. Thus the flow rate can be varied for different test 
conditions. Then the flow is gradually expanded through the tunnel 
cross-section without separation because numerous fine mesh screens are 
encountered by the flow along the way. 

Next, the flow goes through a turbulence management section which 
has two fine-mesh screens, a 12.7 cm length of packed straws, and five 
more fine-mesh screens. When the flow passes through the turbulence sec- 
tion, small eddies are formed which dissipate much quicker than the 
large eddies. The turbulence management section thus keeps the turbu- 
lence level down. 

Having left the turbulence management section, the air enters into 
a contoured nozzle leading to the test section. This axisymmetric noz- 
zle was designed to produce a minimum adverse pressure gradient on the 
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boundry layer to avoid flow unsteadiness associated with local sep- 
aration regions. The area ratio of the cross sections of the turbulence 
management section to that of the nozzle throat Is approximately 22.5. 

The diameter, d, of the nozzle throat Is approximately 15 cm. 

Next, the air enters the test section. The test section Is composed 
of a swirler (optional), an expansion block, and a long plexiglass tube. 
The swirler currently available Is a variable vane-angle type device to 
impart swirl to the flow entering the test section. The expansion 
block, attached after the swirler. Is a 30 cm diameter disk of wcod. At 
present, there are three expansion blocks, and the appropriate choice 
gives a * 90, 70, or 45 degrees. The flow Is expanded into a plexiglass 
tube of diameter, D, of 30 cm, thus giving diameter expansion ratio (0/d) 
of 2. 

A typical real combustor, shown in Figure 1, Appendix B, is ideal- 
ized in the present study, as there are no film cooling holes or dilu- 
tion air holes, and the chamber wall of the test section Is a constant 

diameter pipe. The test section Is carefully aligned using a laser beam 

so that the test section and wind tunnel centerline are colinear. 

2.2 Hot-Wire Instrumentation 

Figure 4, Appendix B, shows the circuit diagram for a constant 
temperature anemometer. The anemometer used for the present study is 
DISA type 55M01 , CTA standard bridge. A normal hot-wire, type 55P01 , 
manufactured by DISA, is used to carry out the measurements of time-mean 
and root-mean-square voltages. These probes have two prongs set approx- 
imately 3 mm apart and carry 5 urn diameter wire which is gold plated 
near the prongs to reduce end effects and strengthen the wire. The 
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time-mean voltage is measured with Hickok Digital Systems, Model DP100, 
integrating voltmeter and the root-mean-square voltage is measured using 
Hewlett Pjckard, Model 400 HR, voltmeter. 

The hot-wire is mounted on the facility with the help of a travers- 
ing mechanism shown in Figure 5, Appendix b. It consists of a base that 
is modified to be mounted on the plexiglass tube of the text section at 
various axial locations. The hot-wire probe is Inserted into the tube 
through a rotary vernier and the base. The rotary vernier is attached 
to a slide which can traverse up to approximately 14.5 cm. Thus it 
becomes possible for the probe to be traversed at any location in the 
combustor fiowfield and rotated through 180 degrees. Figure 6, Appendix 
B, shows the test section with the probe mounted on it. 

L\3 Calibration Nozzle 

The hot-wire is calibrated on a small air jet. The facility con- 
sists of a compressed air line, which delivers the desired flow rate 
through a small pressure regulator and a Fischer and Porter Model 
10A1735A rotameter. The jet housing consists of an effective flow 
management section followed by a contoured nozzle with a 3.5 cm diameter 
throat. 

A rotary table is used to hold the probe while it is being cali- 
brated in three different orientations which are discussed in Chapter 


III. 



CHAPTER III 


STATISTICAL ANALYSIS PROCEDURE 
3.1 Response Equations 

The six-orientation hot-wire technique requires a single, straight, 
hot-wire to be calibrated for three different probe directions in order 
to determine the directional sensitivity of such a probe. The three 
directions and the three calibration curves are shown in Figure 7, 
Appendix B. Each of the three calibration curves is obtained with zero 
velocity in the other two directions. The calibration curves demon- 
strate that the hot-wire is most efficiently cooled when the flow is in 

a 

the v direction. Whereas, the wire is most inefficiently cooled for the 

a 

flow in w direction. Each of the calibration curves follows a second 
order, least square fit, of the form: 

E 2 * A + BZ* + CZ (1) 

where A, B, and C are the calibration constants and Z can take a value 

A A A 

of u, v, and w for the three calibration curves, respectively. 

When the wire is placed in a 3-dimensional flowfield, the effective 
cooling velocity experienced by the hot-wire, in terms of the probe 
coordinator and pitch and yaw factors (6 and K) as defined by Jorgensen 

(20) is: 


Z 2 * v 2 + G 2 u 2 + K 2 w 2 (2a) 



A ^ * t 

u (w, v * 0) 


9 


10 


K „ V (*. u ■ 0) 
w (v, u ■ 0) 


(2c) 


evaluated from the three calibration curves for a constant value of E 2 . 

To carry out measurements in the combustor flowfield, the wire is 
aligned in the flow in such a way that in the first orientation, the 
wire is normal to the flow in the axial direction and the probe coordi- 
nates coincide with the coordinates of the experimental facility. Thus 
the six equations for the Instantaneous cooling velocities at the six 
orientations, as given by King (21) are: 


2 2 2 2 2 2 

Zj 8 v + 6 u + K w (3) 

Z 2 = v ♦ G (u cos 30° ♦ w sin 30°)‘+ k‘ ( w cos 30° - u sin 30°)' (4) 

Z 3 s v ♦ G*(u cos 60° ♦ w sin 60°) + K (w cos 60° - u sin 60°) (5) 

2 2 2 2 2 2 

Za * v + G w + K u (6) 

Z 6 = v + G (w sin 120° + u cos 120°) 2 + if (u sin 120° - w cos 120°) (7) 

Z 6 - v + G‘ (w sin 150° + u cos 150°) 2+ K* (u sin 150° - w cos 150°) (8) 


Replacing the sines and cosines and expanding the square brackets: 


1 

2 

’2 

’3 

$ 

4 

.2 


2 2 2 2 2 

V + G U + K W 

2 3 . W 2 


v" + G (u uw k |) ♦ kVj ♦ J * uw r f) 
v‘ ♦ G : (J-'+ w* l + uw k |) ♦ K 2 (J 2 + u ? J - uw ^) 


2 3 




2 .2 2 22 

V + G W + K U 


2 


V 


2 


V 


+ 


+ 



G (u 


2 3 


w 2 
4" ' 


2, w 2 27 

UW ^) * K (j + U j ♦ UW 

A 2 23 m2 

uw j) ♦ K (w ^ ♦ y ♦ uw 




(3a) 

(4a) 

(5a) 

(6a) 

(7a) 

(8a) 


Solving simultaneously any three adjacent equations provides expressions 
for the instantaneous values of the three velocity components, u, w, and 
v, in terms of the equivalent cooling velocities (Z^ , Z 2 » and Z 3 for 
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example, when the first three equations are chosen). Kin? refers to 
these Instantaneous velocity components as FI, F2, and F3 as follows: 


FI 


F2 


F3 


• 

JaO ♦ (AO' ♦ t') 1 } * 1 

\ T (E^-RT)J 

L x 

| j-AO ♦ (AO 2 ♦ y 2 ) J J 


(6 r - K T ) 


CO - * (AO 2 ♦ 


(9) 

( 10 ) 

(ID 


The values of AO, BO, and CO depend on the set of the thrpe equa- 
tions chosen and are given in Table 1, Appendix A, for appropriate 
equation sets. 

However, these equations cannot be directly used because it is im- 
possible to obtain Zj, Z^, and Z ? at a single instance in time. There 
fore Equation 9 through 11 must be expressed In terms of mean and root- 
mean-square values. Equation 1 can be written as: 

<t>( E i ) - Z i - |j-B ♦ Jb* - 4C (A - E^J* ]/2cj * (12) 

The above equation is in terms of Instantaneous velocity Z^ and instan- 
taneous voltage E^. In order to obtain an expression for time-mean 
velocity as a function of time-mean voltage, a Taylor series expansion 
of Equation 12 can be carried out. 


Since Z i ■ *(E,, + Ej) ( 2 

l \ a + e!) = ♦(ty + jf- 


( 13 ) 


ORIGINAL P"'*: r 12 

OF POOK QUALITY 

The Taylor series Is truncated after second order terms assuming the 
higher order terns to be relatively small. Time averaging both sides of 
the above equation and employing the fa'.i. that £'■ 0, yields: 



where $ Indicates that the function Is evaluated for E^. To obtain 

2 2 

IJ * o y , the relationship as given by Hlnze (13) Is: 

1 4 1 

Zj 2 « Oj 1 * Expec [Z 1 2 ] - (Expec [Z^]) 2 

— 2 

Since Expec [Z|] - <t> ♦ 1/2 -y 2 * °E i ’ 

the differential in Equation 16 can be evaluated as: 

!?f ■*<#,> 2 

Then Equation 16 becomes: 

Expec [Z-.|] “ e 2 ♦ (^- ) 2 • • o| 

Squaring Equation 14 and substituting with Equation 18 into Equation 15 

gives: 

■ v • f 2 - "s,* - <'« ^ • ° E ;>’ as) 

Thus Equations 14 and 19 give the mean and variance of individual 
cooling velocities in terms of the mean and variance of the appropriate 
voltage. 

In a 3-dimensional flow, It Is usually desired to obtain the mean 
and variance for the individual velocity components in axial, azimu- 
thal, and radial directions, and also their cross correlations. 

The procedure to obtain the mean and variance of the individual 


05 ) 

(16) 

(17) 

(18) 
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velocity components Is the same as for the effective cooling velocities 
except that u, w, and v are functions of three random variables and 
there are extra terms In the Taylor expansion to account for the co- 
variances of the cooling velocities. Thus the three mean volicltles as 
given by Dvorak and Syred (19) and King (21) are: 


i j • 3*F1 

“ ' fl(z P- z q- z R' * 7 3T757J • 


( 20 ) 


where time-mean values are to be understood on the right side of this 
and subsequent equations. 


3 a F2 


" * f 2U p .Zq.Z r ) * ? * } 57^ ' °Z 


* J F 2 


3 
1<J 


Vj 


( 21 ) 


and 


- n /1 7 7 » x 1 J 3 2 F3 
v - F3(Z p .Z Q .Z R ) ♦ 7 E jy-T 


1-1 1 


3 

Z 

1<J 


3 : F3 


Vj 


( 22 ) 


where K ? 7 is the covariance of the cooling velocity fluctuations and 
' j 

Is defined as: 



(Zj - r,)(Zj - ?,>« 


(23) 


Also the normal stresses are given as: 
-tt l 1 3FK t 2 * I l 3H 3F1 

1-1 3Z 1 z i 1 j 3£ i dL i 

IFj ^ 

r 3 ' n . k 
1<j * <Z 1 Z j 


-E (f^)’ 

1-1 3Z 1 


3 3 -F? 

0 , 1 ♦ I E 

JLj 


■1 


3 2 F2 


1 j *“1 

W - 


3F2 

57J 


t<“j ^ ’ z l z j 



1 t. 3*F2 

7 1 *7^ 
1 1-1 




(24) 


(20) 


2 
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and 


■ i &v 

1-1 3Z 1 


. . )F3 3F3 _ 

. oZ< ♦ E £ 57“ • 377 

1 i j 


l 3 2 F3 
1<j 


ir«J 


Vj 


1 l 3 2 F3 
1-1 ot 1 


7 ,5, 577 • °z, * 


(26) 


Also the shear stresses as given by Dvorak and Syrcd (19) arc: 


(TV 


3 3 


♦ EE 


i j 
1<j 


* L 


3 3 
E E 


1<j 

and finally. 


: 3Fi 

3F2 

2 . 1 l 3F1 

3F2 , 


1 

i ^ • 

327 • 0 

■*1 

♦ EE . 7 — . 
1 j 3Z 1 

^ ‘ Z 1 Z j 


7 



W 



— 

3* FI 



1 I 3 a F2 

3 3 

_ 2 i r r 

3 2 

F2 

57752J 

Vj 


? 1-1 * 

0 j ^ E 0 

z i 1 j 

52737 , 




- 

i<j 



: 3Fi 

3F3 

2 ♦ r ? 3F1 

3F3 „ 


1 

1 3Z i * 

3T, • c 

■*1 

♦ EE xy . 
1 J 3Z 1 

377 • Vj 


7 



1 |*j 



mm 

3 2 F1 



■» 

1 r 3 2 F3 

3 3 

_ 2 » r r 

3 2 F3 

W] 

* K Z Z 
z 1 z j 


? ,£, 577 • 

Oy ^ L. Y. 

Z 1 i j 

¥J 




»• 

1<J 




l 3 * F1 a 

i«i 3Z i L \ 


Vj 


(27) 


1 * 3 2 F1 _ 2 

1 W 7 ‘ °7 
1-l 3Z i Z 1 


%l. 


(28) 


-r-r l 3F 3 

w 7 v - £ 57 " • 377 

i*l az 1 1 


33 3 2 F2 


3 3 aF2 3F3 

% + f J 3^ ' ^ 


j Z i 


L L 57“a7~ ’ Z,Z- 
i j 3Z 1 3Z j 4 1 4 j 

i<j 


1l*j 
1 1 li£3 


1 J 3 2 F2 2 + 

* ,£,*7 ' z i 


a 2 £3_ 

+ L Z 3Z ; 3Z ; 


■1 


1 j i J 
1*j 


Vj 


( 29 ) 
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3.2 Calculation of Covariances 


Dvorak and Syred (19) used a DISA time correlator (55A06) to find 
the correlation coefficients between the velocity fluctuations In tiie 
three directions. The method adopted by Klnq (21) Is to use the Infor 
matlon obtained by all six orientations and devise a mathematical pro- 
cedure to calculate the covariances. 

The covariance matrix as derived by King Is: 



where 



(30) 


(31) 


n 


2 



1 

7 


P .2 

- 3p 


1-p 


3Z 


2 

7 



2 


(32) 
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and 


R a Z 

1 E ^ °3 

"3 ' ! P*i • °3 * i t-p • 0 z , 2 


Al SC 


7 2 ,7 ^ ► 27 2 

*P “r 


p 2 • 22/ ♦ *„* * Y 


2 ♦ 7 2 

3 ~ “P ' W -Q R 


and p, * 2lJ - 2l n c ♦ Z 


(33) 


(34) 

(35) 

(36) 


Matrix (T) Is a three by three matrix and is given In Table II. 

King discovered that matrix (T) Is a singular matrix for all cases 
and hence equation 30 cannot be solved. Therefore, to get covariances 
one needs extra information. King has made an assumption about the 
relationship between the covariances In the form: 


^ z * *z z 

W n -Y 11 R 


(37) 


rfhere n is given a numerical value of 0 . 8 . 
Also K_ 


ZpZg s obtained from the quadratic equation: 


Vq 


b27 p ? -n 


t Kz ^ 

* Vq 


6 I n T n - 


P Q l 




(ni-Z; 


3 

R+l 


•‘ 3 * Z R+3 ) 


♦ fni« 


^R+l " 2 n 3 *^R+3 


» 0 


(38) 
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Equation 38 provides tho two values for K . The covariance is 

Z P Z Q 

related to tl.e correlation coefficient as: 


Z P Z q 


Vg [V .Vi 1 


(39) 


whe»e -1 < y * 1 

z r j 


Therefore, tquation 39 is written in the form: 



The two calculated values of K from equation 38 are then subst i ited 

Z P Z Q 

in Equation 39, and the two correspondinq values of y are calculated. 

Z P Z Q 

The correlation coefficient which lies within the required range of 
+ 1,1s used. For the case when the absolute values of both the correla- 
tion coefficients are larger than 1, the covarijnce is given by 


C Z Z * 0,9 [V * °z 2 1 * 
Z P Z Q l P qj 


(41) 


Having calculated K , K can be calculated from the relationship: 

Z P Z Q Z Q R 


K, 




z„z„ ' rxx: I 2 - Wz p z 0 * "i- z p,3 - "3 


■ 

s 

• Z p*5 


(42) 


1 I _ «■ _ _J 

= y 

Q'-R ‘■Q t R 

A similar test is applied to ensure that the absolute value of y is 

Z Q Z R 

less than one otherwise K is calculated from the relationship: 


Z Q Z R 


%Z R ■ °- 9 [v • v] 


(43) 
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K can now be calculated from equation 37. The calculated value of 

Vr 

K , K , arid K can now be substituted in equations 20 thru 22, 
Z P Z Q z Q 7 R z P z R 

and 24 thru 29 to calculate the mean velocities and Reynold stresses. 

It was observed during the present study tnat King's method Is not 
self-consistent In calculating the covariances. The correlation co- 
efficients were found to have values greater than one and therefore It 
was necessary to have a more consistent method to calculate the covari- 
ances. Occasionally, King's method assumed that y and y had 

Z P Z Q Z Q Z R 

values of 0.9 and y had a value of 0.648. But this was done only 

Z P Z P 

when some of the correlation coefficients were greater than one. The 
present method assumes constant values of the correlation coefficients. 
King has suggested that If two wires are separated by an angle of 30 
degrees, the fluctuating sianals from the wires at the two locations 
would be such that their contribution to the cooling of the wire would be 


related by the cosine of the angle between the wires therefore, y 
cos 30° * 0.9 and similarly we would get 



y Vr " °- 9> 

alS0 1z p z R ■ n* Vo * % Z R • °- 648 

Therefore the present method allows the covariances to be calculated 
using the following three equations: 



’HAPTER IV 


UNCERTAINTY ANALYSIS 

An uncertainty analysis Is presented here with a view to demon- 
strate the reliability of the six-orientation hot-wire technique and its 
sensitivity to various input parameters which have major contributions 
in the response equations. The analysis Is done for both laminar and 
turbulent flow cases. The salient results are tabulated in Tables III 
and IV of Appendix A. 

4.1 Effect of Pitch and Yaw Factors 

Pitch and yaw factors (G and K) are used in the response equations 
described In Chapter III in order to compensate and account for the 
directional sensitivity of the single hot-wire probe. Figure 8, Appen- 
dix B, shows the pitch and yaw factors plotted against the hot-wire 
mean elective voltage. Both the pitch and yaw factors are functions of 
the hot-wire mean effective voltage, but the yaw factor is far more sen- 
sitive. A 10 percent increase in the voltage reduces the yaw factor by 
56 percent and the pitch factor by 13 percent. The value of the pitch 
factor stays very close to one and hence does not have a major contribu- 
tion in the response equations. For this reason, it is necessary to 
further consider the yaw factor, which is now examined for both laminar 
and turbulent flow conditions. 
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4.1.1 Lam inar _ Flow 

For laminar flow cases, the covariances K become zero and d op 

Z i Z J 

out of the response equations. Then equations 20, 21, and 22 can be 
written as: 

u = F1 ^p’ *Q* Z r) 

w • n <V 2 q- z r> 

» ■ F3 < Z P - V z r.> 

Experiments were performed on a calibration nozzle free jet in the po- 
tential core where the flow can be idealized as being laminar. 

As Table 111, Appendix A, shows, the effect of yaw factor on time- 
mean axial and swirl velocities is insignificant for the laminar flow 
case. 

4 . K 2 Turbulent F low 

The variation of yaw factor is studied on the turbulence quantities 
such as mean velocities, turbulence intensities and the shear stress 
u'v'. As stated in Table III, all turbulence quantities behave differ- 
ently to the variations in the yaw factor. The effect on all the turbu- 
lence quantities, except the mean radial velocity, is insignificant. In 
the case of mean radial velocity* the term (G -K ) in the denominator of 
Equation 9 changes the value of F3 considerably for small changes in the 
yaw factor. 


4.2 Effect of Correlation Coefficients 


Correlation coefficients are used in Equation 40 to calculate the 
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covariances between the fluctuations of the cooling velocities experi- 
enced by the hot-wire at adjacent orientations. These are then used In 
Equations 20 through 29 to calculate various turbulence quantities. A 

wide range of correlation coefficients (>_ , ) between 0.1 to 0.9 are 

Z P Z Q 

used to study the behavior of the turbulence quantities. Among nil ihe 

turbulence quantities, u'v’ was found to be most sensitive to variations 

in the correlation coefficient (v ). In view of the sensitivity of 

Z P Z Q 

u'v' to > and the assumptions required to estimate > , it is 

Z P Z Q Z P Z Q 

apparent that this is the major source of the significantly large un- 
certainty in the estimate of the turbulent shear stress. This appears 
to be an inherent deficiency of the six-orientation single hot-wire 
method. 

King (21) used a parameter Eta (n) to relate the covariances be- 
tween the fluctuations of the effective cooling velocities that are sep- 
arated by 30 degrees with the covariance of velocities separated by 60 
degrees (see Equation 37) He suggested a numerical value of 0.8 for n. 
Table III shows the effect of n on the turbulence quantities to be in- 
significant and hence the present study retains this value of 0.8 in all 
subsequent deductions. 

4.3 Experimental Uncertainty 

Experimental uncertainty was tested for both laminar and turbulent 
flow cases. The main reason for these tests was to determine the mean 
and variance of the output quantities when obtained from the six possi- 
ble choices of three from among the six possible response equations 
(Equations 3 through 8 in Chapter III). Another objective of the study 
was to judge the extent of errors in output quantities because of errors 
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in measurement of mean and root-mean-square voltaqes. 

4 . 3 J L aminar FI o w 

The calibration free jet facility was used to conduct laminar flow 
uncertainty experiments. To generate velocities In the axial and azimu- 
thal direction with respect to the wire, the wire was offset by 45 de- 
grees to the main direction of the flow and placed in the potential core 
region, thereby achieving two equal components of axial and swirl velo- 
cities. However, upon data reduction. It was observed that the two com- 
ponents were not equal. The variation among the two components was 
different for each choice of the six combinations of three ddjace.it 
response equations. In general, the variation among the two components 
was negligible. 

Table III shows the effect of variations in measurements of the 
hot-wire mean voltages on the turbulence quantities. For laminar flow 
case, the mean axial and swirl velocities are extremely sensitive to 
errors in measurements of hot-wire mean voltage. This particular test 
stresses the need for using precise voltmeters. A 10 percent error in 
measurement of one of the six mean voltages leads to an error of 90 
percent in axial velocity deduction for the conditions of this test. 

At other flow conditions, similar gross sensitivity may be expected. 

Turbulence quantities are calculated using six different combina- 
tions ot the three mean effective cooling velocities experienced by 
the hot-wire at three adjacent orientations. Table IV, Appendix A, 
demonstrates good consistency between the six possibilities for mean 
axial and swirl velocity determination in laminar flow conditions. 



23 


4 . 3 J 2 T ur b ulen t _F 1 ow 

As observed for the laminar flow case, errors in mean voltage 
measurements are extremely magnified in calculations of turbulence 
guantities. Table III shows these large variations in the turbulence 
quantities. 

For turbulent flows, a large scatter Is observed among the six 
values of turbulence quantities deduced from the six different combina- 
tions. To get an estimate of the scatter, the flowf.eld location x/l) 
2.0, r/D * C.25 for the case of side-wall angle a * 45 degrees was 
selected inside the main test facility. At this location in the flow- 
field, the turbulence quantities obtained are good representatives of 
turbulence level in the combustor flowfield. 

Table IV shows that for turbulent flow, all the six combinations 
do not reveal all the turbulence quantities. The omitted i vems 
correspond to occasions when the velocity function f'3 attains a complex 
value via the requirement of the square root of a negative value. Then, 
no further prooress could be made with that particular sec of three 
adjacent orientations in such situations. 

Table IV also highlights the scatter among the six values of each 
turbulent quantity when solved using six different combinations. It is 
evident that certain quantities (such as mean radial velocity, the 
radial turbulence intensity, and the shear stress u'v') have very large 
scatter. This shov/s a great uncertainty in the use of six-orientation 
hot-wire technique in measurement of these quantities. 


CHAPTER V 


RESULTS 

The six-orientation hot-wire technique is employed to measure the 
turbulence quantities for nonswirling conditions. The experiments have 
been conducted for expansion angles of 90 degrees (sudden expansion) and 
45 degrees (gradual expansion). A computer program, listed in Appendix 
D, written in Fortran language, is used to process the data on an IBM 
370/168 computer. For each location in the flowfield, six combinations 
of three adjacent orientations are selected and six values of each of 
the nine turbulence quantities are obtained. So, a decision has to be 
made about the selection of one of the six values. In nonswirling con- 
ditions, the flow is mainly dominated by the axial velocity. When the 
hot-wire is parallel to the axial direction; it experiences the least 
cooling effect from the axial velocity, whereas the radial and swirl 
velocities affect the wire most efficiently. Therefore, a small change 
in the v and w velocities will show a significant change in hot-wire 
voltage. Hence the set of orientations labeled (3, 4, 5) in Chapter III 
(orientation 4 having the hot-wire parallel to the x-direction) is 
chc en and used in all subsequent results presented, except whtre noted 

otherwise. Nevertheless, there are some quantities, such as v' and 

rms 

u'v', which appear to be better represented by alternative sets of three 
adjacent orientations, but the appropriate choice is not known a priori. 
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5.1 Mean Velocities 

Radial distributions of time-mean axial and radial velocities <»re 
plotted in Figures 9 and 10, Appendix B, respectively. Mean axial 
velocities for different axial locations and expansion angles are com- 
pared with measurements done with a crossed hot-wire probe by Chaturvedi 
(5;. A good agreement is found between the two studies. 

Because of the inabil :y of the hot-wire to determine the sense of 
the flow direction, the presence of the corner recirculation zone was 
observed by a sudden increase in the axial velocity closer to the wall. 
Mean radial velocity was found to increase at the centerline with an 
increase in the axial distance. The mean velocity profiles tend to 
get flatter further downstream from the inlet. For a « 45°, mean radial 
velocity at the centerline increased from 5 percent of the maximum in- 
let mean velocity at x/D * 0.5 to 16 percent of the maximum inlet velo- 
city at x/D ■ 2.0. A similar increase was observed for c* = 90°. 

5.2 Turbulence Intensities 

The six-orientation hot-wire technique enables one to measure the 
axial, radial, and azimuthal turbulence intensities at various axial 
and radial locations in the confined jet flowfield. The radial distri- 
butions of these turbulence intensities are plotted in Figures 11, 12, 
and 13 of Appendix B. The axial and radial turbulence intensities are 
compared with Chaturvedi 's study (5) and reasonable agreement is found 
in the case of axial turbulence intensities. However, the two studies 
are not in good agreement for radial turbulence intensities. The peak 
values measured in the present study are much lower, in certain cases 
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being only 50 percent of the previously measured peak values (5). While 

solving the r ix sets of combinations of three adjacent orientations, it 

was found that v’ /*j has a large scatter. For example at x/D ■ 2.0, 
rms o 

r/D 0.300, and x ■ 45°, the mean and the standard deviation of v' / 
u () , among the six sets of readings, were found to be 0.1447 and 0.0330, 
respectively. This large scatter shows that In nonswirling flow this 
technique is not a very accurate way of measuring the radial turbulence 
intensities. Nevertheless, results shown in Figure 11 have been obtain- 
ed with the set of orientations (3, 4, 5) being used. 

5.3 Shear Stresses 

In nonswirlinq flow conditions, measurements have been made of the 
turbulent shear stress u’v'. The radial distribution of u’v'/u^ at 
various axial locations is plotted In Figure 14, Appendix B, and is com- 
pared with the earlier study done by Chaturvedi (5). In general, the 
two studies are in good agreement, but they do differ in two respects: 
the centerline values far downstream and the peak values near the inlet. 

Chaturvedi (5) measured u^v' to be zero at the centerline at all 
axial locations. However, in the present study, u'v' is found to be non- 
zero at the centerline at axial locations greater than x/D r 0.5 for 

both side-wall angles o * 45° and a = 90°. A detailed study shows that 

? 

the scatter for u v /u , when calculated from different sets of adja- 
cent orientations, is quite large. The ratio of standard deviation to 
the mean is approximately 0.6 and varies with position. 

Peak values of u'v' are seen to be in good agreement except close 
to the inlet. At x/D = 0.5, Chaturvedi (5) measured peak values 
approximately 50 percent higher than in the present study. It must be 
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remembered that there is always difficulty In measuring shear stress 
values in thin shear layer reqions. In the present study, there is 
also the previously-discussed deficiency, see Chapter IV, because of 
the assumpitions made about the correlation coefficients 7 

Z 1 Z j 

These assumptions may be the major source of significantly large un- 
certainty in the calculation of turbulent shear stress values. 



CHAPTER VI 


CLOSURE 
6.1 Summary 

The six-orientation hot-wire technique Is a relatively new method 
to measure time-mean values and turbulence quantities In complex three- 
dimensional flowfields. Applied In this study to nonreacting nonswirl- 
ing axisymmetrlc flowfields, measurements of time-mean and root-mean- 
square voltages at six different orientations contain enough Information 
to obtain the time-mean velocities, turbulence intensities and shear 
stresses. At each location in the flow, there are six different values 
of each of the above quantities that can be obtained by using six sets 
of measurements of three adjacent orientations. Because of axial velo- 
city domination, a particular set of orientations was chosen. Neverthe- 
less, the measurement accuracy can be well judged by the scatter of the 
values of turbulence quantities among the six different combinations of 
sets of three mean effective coolinq velocities. The nonswirling 
confined jet flow was investigated with this technique. It was found 
to be an excellent method to find time-mean velocities. It also gave 
good results for turbulence intensities and shear stresses. An uncer- 
tainty analysis done on this technique reveals that certain output 
parameters such as the axial, radial, and azimuthal turbulence inten- 
sities and shear stresses are extremely sensitive to some input 
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parameters such as yaw factor and mean voltages. 

6.2 Further Work 

The multi -orientation sinole-wire technique is a useful cost- 
effective tool for the investigation of complex flowfields. At present, 
there is a need to check repeatability under nonswirling conditions 
before progressing to the investigation of flows with moderate and 
strong swirl. This would lead to further evaluation of reliability and 
accuracy of the technique in general flowfields. Thus far, there is an j 
priori assumption about the evaluation of covariances, which entails the 
use of constant values for the correlation coefficients. Further work 
might call for the development of alternative methods to specify the 
covariances. Nevertheless, the method has potential for further use in 
the experimental evaluation of complex flowfields. 
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TABLE II 

Matrix (T) in Equation 30 


2 Z P Z g Z P Z R 4 Z Q Z R 

Z P+3 Z P+3 Z P+3 
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4 i p i R 


P+4 


?Z Q Z K 
Z P+ 4 


4 Z P 
Z P+5 


2 Z P Z R 
Z P+5 


2Z Q Z R 
Z P+-5 



TABLE III 
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. Typical Axisymmetric Combustion Chamber of a 
Gas Turbine Engine 


Figure 1 





(b) EXPECT ED KECIRCWLATION ZONES 


Figure 2 . The Flowfield Being Investigated 
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Figure 3. Schematic of Overall Facility 
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Figure 5. Manual Traversing Mechan- 
ism Used for Hot-Wire 
Orientations in the 
Flowfield 




Figure 7. The Three-Directiona - * hot-wire Calibiatior 




PITCH , YAW FACTORS 
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HOT-WIRE VOLTAGE E(VOLTS) 


Figure 8. Plot of Pitch and Yaw Factors Versus Hot-Wire Voltage 







Figure 12. Radial Distribution of Radial Turbulence Intensity in Nonswirling Confined Jet 






Figure 13. Radial Distribution of Azimuthal Turbulence Intensity in Nonswirling Confined Jet 
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USf R 'S GUIDE TO COMPUTER CODE FOR SIX-ORIENTATION 
HOT-WIRE DATA REDUCTION TFCHfllOUE 

A canputrr code i< developed to obtain the turbulence quantities 
* the te ique discussed In Chapter III. Measurements in a turbu- 
i 'wf’eld orain six mean and six root-mean-square voltages. A 

thr- 11 recti one 'Ot-wire calibration reveals three calibration con- 
stan n each d ction. The input to the computer code Is the mean, 

no ro-nt-nx ,n-squan vc • iges and also the calibration constants. The 
-erimt tal data then processed by the MAIN subprogram and various 
; '-out ini to get ie output in the form of nine turbulence quantities 
ting 'f the three mean velocities, the three turbulence intensi- 
and th three shear stresses. To facilitate the use of the com- 
ite code, tl function of each subprogram is discussed here in detail. 

he MAIN Sudi rogram 

MAIN is the major part of the computer code which accepts the input 
in the form of mean and root-mean-square voltages and calibration con- 
stants and calls various subroutines to solve the equations listed in 
chapter III and finally calculates the turbulence quantities. 

( 8 ) Calculation of Mean Effective Cooling Velocities and 
Variances 

Main calculates the six mean effective cooling velocities 
using Equation 14. This equation employs the input values 
of six mean voltages and calibration constants in u- 
direction (see Figure 7, Appendix A). The MAIN then cal- 
culates the six values of variances using Equation 19. 
Equations 14 and 19 give mean and variance of individual 
cooling velocities in terms of the mean and 



variance of the appropriate volt aye. 

Calculation of V elocity Function s and Differ enti a is 
Having calculated the mean effective cooling velocities 
and variances, the MAIN then calls various subroutines 
to obtain the necessary Information required to calcu- 
late velocity functions using Equations 9 through 11. 

The main then calculates the first and the second dif- 
ferentials of the three velocity functions with respect 
to the three s* t itean effective cooling velocities. 
The different . are given as: 
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(111) Calculation of Covaria nces 

At this stage, the user has the option, wheth< 
dilate the covariances oy using King's (21) method or 
by assuming constant values of correlation coefficients. 

To get the covariances using King c method, the MAIN 
has to call the subroutine COVAR, otherwise MAIN calcu 
lates covariances using Equations 35, 41, ana 44. 

( 1 v ) C alculation of the Turbulence Quantitie s 

Now the MAIN has all the information needed to calculate 
the mean velocities using Equations 20 through 22, also 
to calculate the tu. oulence intensities usino ‘•-.ur lor 
24 through 26, and finally to calculate the shear 
stresses using Equations 27 through 29. The MAIN then 
prints out the normalized values of the turbulence 
quantities in the form of nine two by three matrices each 
containing the six values of a turbulence quantity cal- 
culated using six different combinations. 

2. Subroutine CPYF 

This subroutine calculates the pitch and yaw factors using the cal- 
ibration constants obtained by three-dimensional calibration. Tht equa- 
tions used to calculate these factors are: 

e ■ i Si 

u (w, v * 0) 

I/ - V (w, U - 0) 

K w (v, u - 0) 

evaluated at a constant value of E 2 . u, v, and w are obtained using 
equation 12 for their respective calibration constants. The value of E 
can be adjusted to obtain an interval aE to get appropriate values of 
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G, and K. 

3 . Subroutine FMCV 

The task of this subroutine Is to find the mean effective oling 
velocity which has minimum value among the six calculated b\ t MAIN. 
FMCV also finds the two mean effective cooling /elocit1< > which »re 
adjacent to the minimum mean effective cooling velocir nd returns the 
set of the three to be used by MAIN for further data pn cessing. 

4. Subroutine SEA8C 

SEABC recognizes the three selected mean effective cooling v< ji i- 
ties Zp, Zg, and Z R , and sets the three appropriate equations for AO, 
BO, and CO In terms of Z p , Zg, a. id Z R , using Table V. AO. “0. «nd CO, 
are used by M. N to calculate the three velocity functions given t 
Equations 9 through 11. 

5. Subroutine CDABC 

CDABC calculates the first and second differentials of AO, BO, and 
CO with respect to Z p , Zg, and Z R . It is evident from Table V that AO, 
BO, and CO are functions of Z p , Zg, and Z R and so are their first and 
second differentials. 

6. Subroutine COVAR 

This subroutine calculates covariances using a method suggested by 
King (21). This method calls for employing Equations 40 through 43. 
This subroutine can be called only when one desires to calculate covar- 
iances using King's method. Otherwise, the covariances are calculated 
within MAIN by the procedure already described. 
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TABLE V 

LIST OF FORTRAN VARIABLES AND THE 1R 
MEANING IN RESPONSE F QUA! IONS 

Input Value? 

EM Z 

ER E'rms 

Mean Fffpctive Cooling Velocities and Variam.es 
AMECV l 

VAR 0 2 J 

Pitch and Yaw factors 
PF G 


YF 


K 


Set of Three Cooling Velocities Chosen 


IP 

ZQ 

ZR 


Derivatives of Functions AO, BO. 


A1 

|A0 


3Z P 

A 2 

MO 


3Z o 

A3 

?A0 



A21 

a 2 AO 


V 

A22 

a 2 ao 


3Z q 2 

A23 

a 2 AO 




and CO (Definitions for B and C 
analogous to tl 
defined here) 


are 
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TABLE V (Continued) 
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iG 


az p 

DF1Q 

3_f7 


3Zq 

DF1R 

iZI 


3Z r 

D2F1P 

a ? rr 


aZp 2 
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3 V z i 


Covariances 


TABLE V (Continued) 


Output . .11 iables CaU i1at< d 
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LISTING OF THE COMPUTER PROGRAM 
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ABSTRACT 


The general aim of the present study is to investigate ax' symmetric, 
nonreacting flow in sudden expansions. This configuration simulates the 
typical geometries found in can-type * 2 ? turbine and ramjet combustion 
chambers. The turbulence quantities reported here are used to understand 
the complex flowfield in question. They are also to be used to recomnend 
appropriate turbulence model advances in the simulation of combustor chamber 
flowfields. In the present contribution, data are obtained with a hot-wire 
probe, with two inclined wires, and specialized electronics capable of 
instantaneously adding and subtracting the two voltages. 

Initially measurements were taken in a free shear flow, in the fonn of 
a free ax i symmetric jet. The main results reported are properties of the 
axisymmetric confined turbulent jet in the sudden enlargement [side-wall angle 
or 90 deg] test facility. Experimental measurements of mean and fluctuating 
velocities are presented at several axial stations under nonswirling condi- 
tions. Results are compared with those of previous researchers where possible, 
and good agreement is to be noted. 
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1. INI RODUCTION 


A major objective of the combustion aerodynamics research at Oklahoma 
State University Is to perform detailed turbulence measurements In an iso- 
thermal swirling confined jet. The geometry is a simplified simulation of 
an aircraft can combustor. This report discusses measurements made with 
a crossed hot-wire probe in this flowfield. 

Two independent hot-wires on one probe Inclined at ♦ 45° to the f.ow 
direction were used as early as 1951 by Corrsin and Uberoi,* where they 
measured the power spectra of velocity and temperature fluctuations in 
heated and unheated jets, further work by Davies et al ? and Koplin 3 con- 
sisted of using single and double crossed wires in the measurement of 
turbulence intensity and shear stress in tne mixing region of a round jet. 
Corrsin and Uberol 1 *'* and Wygnanski and Fiedler 5 performed measurenents in 
the self-preserving region of the jet. More recently, Bradshaw et al‘ 
measured turbulence properties in the noise producing region of the jet, 
that is, up to 7.5 downstream of the jet exit. As all jets have somewhat 
different nozzle configurations, this study concentrated on free jet 
measurements tar downstream of the jet exit in the self-preserving region 
of the jet. The data of Corrsin 7 were used as a basis of comparison of 
the present hot-wire measurements to provide an independent check on the 
data acquisition and interpretation technique. Measurements of mean and 
fluctuating properties were taken at an axial distance of 20 diameters 
downstream of the jet exit. 

Several research groups have made measurements in axisymmetric , sudden 
expansion flowfields similar to the one being studied.' 10 Ha Minh and 
Chassaing 8 made pitot probe and hot-wire measurements in various sudden 
expansion flowfields. The expansion ratios D/d, they used were 0.6, 2 and 


1 


1.6 out the axial locations they measured at were not compatable to those of 
the present study. Ha Minh and Chassaing* successfully managed to measure 
the Reynolds stresses using the rotating, inclined, single hot-wire techni- 
que. Moon and Rudinger' measured mean velocities and recirculation region 
geometries with the use of Laser Doppler anemonietry equipment but the ex- 
pansion ratio differed from the one used in the present experiments. 

Chaturvedi 10 performed mean velocity and turbulence measurements in an axi- 
synmetric suoden expansion flowfield with an expansion ratio D/d r 2, which 
is identical to that used in our experiments. Consequently his measurements 
can be compared directly with those made in the present study. Also comparable 
are single-wire six-orientation measurements, 11 recently obtained in the same 
test facility as used in the present study. 

Experimental facilities and instrumentation are described in Section 2. 

The measurement procedure is dealt with in Section 3, which includes calibration 
and data reduction techniques, and preliminary proof of principle measure- 
ments i r. a free jet with results compared with Corrsin.' Confined jet 
measurements are given in Section 4, where validation of techniques is further 
demonstrated by way of comparison with earlier results of Chaturvedi 1 ' and 
Janjua. 11 
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? . E XPERI MENTAL FAC1L1T1 E_S 
2.1 ConfJ ned Jet Flow f leld 

The experiments have been conducted In the confined jet facility 
shown schematically In Fig. 1. The faclltly has an axial flow fan whose 
speed can te changed by altering the drive pulley combination. Numerous 
fine screens and straws produce flow In the settling chamber of low tur- 
bulence Intensity. The contraction section leading to the test section 
has bt?n designed by the method of Morel 1 * to produce a minimum adverse 
pressure gradient on the boundary layer and thus avoid unsteady problems 
associated with local separation regions. The sudden expansion consists 
of a 15 cm diameter nozzle, exiting abruptly Into a 30 cm diameter test 
section. The substantial size of this test model will provide good probe 
resolution for the hot-wire measurements. The test section 1r constructed 
of plexiglass to facilitate fiow visualization and ease of location of 
measuring probes. 

Flow at abrupt expansions is a typical example of the problem of 
separation of a flowing fluid from the boundary. The flow Is separated 
at the sudden expansion and causes a region of reverse flow. This is called 
the corner recirculation zone (CRZ). This region is associated with a large 
pressure drop which adversely effects the performance of a can combustor. 

On the shear layer between the two regions the turbulence levels are high, 
which would indicate that good combustion could te achieved at this position. 
Reacting flows in actual can combustors would have points of high heat trans- 
fer at the area of impingement of the dividing streamline. This, of course, 
could lead to rapid deterioration of the combustor and is therefore of great 
interest to combustor designers. 
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2.2 free Jet Calibration Facility 


It was necessary to calibrate the crossed hot-wires In a flow of 
known characteristics. This was achieved by means of a round free jet 
Issuing horizontally Into a quiescent, atmosphere. The jet wis formed by 
using a seamless contoured nozzle of 34 rim exit diameter, d The nozzle 
was fed by a thermally stablllzel compressed air generator which was 
regulated by an upstream control valve to provide the correct flow rate 
through the calibration nozzle. The sensitivity to angulation was oo- 
talned by usage of a rotary table (see section 2.2) on the jet centerline 
to rotate the probe from -10° to +10° to the free stream flow direction. 

The usual hot-wire instrumentation was used to calibrate the probe as 
described in section 1.3. 

2 . 3 Hojt-Wi re Instrumentation 

The probe used in tnis study consisted of two independent inclined 
hot-wires. The wires were inclined at ♦ 45° to the flow direction and were 
made from 5 pm tungsten wire soldered onto the tips of round jewellers 
broaches. The t*. oaches were supported by a pe-forated ceramic insulator. 

The mean and fluctuating flow measurements were made using separate 
constant temperature anemometer electronics for each of the hot-wires. The 
anemometers consisted of a DISA 55M01 main frame with a 5SM10 standard 
bridge. In addition to this it was necessary to instantaneously add and 
subtract the wire rms outputs This was achieved by the construction of 
commercially available integrated circuits. The multiplication of these 
voltages was obtained by using a Saicor model SAI 43 correlation and 
probability analyzer. The fluctuating signals were also amplified, with 
a Hewlett Packard amplifier with 20 dB gain on each of the signals. Normally 
the frequency response of the hot-wires and the associated electronics was 
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approximately 40 kHz based on square wave response test*. A schematic of 
the electronics layout Is shown in figure 2. 

Experience In the use of crossed hot-wire has been gained at Oklahoma 
State University by researchers such as Swearingen** and Morrison. 1 * 
though the majority of their measurement dealt with compressible flow, 
knowledge of the crossed-wire technique was first developed in subsonic 
jets. It was through their previous experience that the use of llnearizers 
was judged unnecessary for the measured quantities in the report. 
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3. EXPERIMENTAL TECHNIQUES 
3 . 1 Data Reluct Ion 

Following the Ideas of Corrsln and Uberol 1 the Instantaneous velocity 
fluctuations measured from an inclined hot-wire can be represented by the 
following expression. 


e' 
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A_ H- ♦ A.. *- 


m u — v u 

where and A y are the sensitivity coefficients for axial velocity tluctua- 
tions and radial velocity fluctuations respectively. The coefficient*. c u n 
be evaluated as follows: 
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The two oppositely inclined wires of a crossed hot-wire probe give 
cutputs that can be instantaneously added and subtracted to result In: 
(el ♦ el) ■ (E.A^ ♦ £,%,,) * (E,A y) - E,A y; ) J' 

(el - ei) - (E,A mi - E.A^) u -' * (E.A V> ♦ E,A yj ) jf 

The crossed hot-wires ate suitably matched so that the sensitivity 
coefficients A (j) and A v are approximately the same for each wire. These 
simplifications result in two voltages proportional to u'/u and v'/u 
respectively 


where 


and 
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where 

»♦ ■ (ElA v, 4 E « A „> 

Additionally, the major shear component of the Reynolds stress tensor 
In a shear flow may be determined by the multiplication of the voltages 
proportional to u'/u and v'/u and using the following relation 


3 . 2 Cal ibration 

The calibration Is performed by placing the crossed hot-wire piobe on 
the free jet calibration facility centerline near the nozzle exit within 
the potential core region. The axial velocity Is then varied by an upstream 
control value and a Rotometer Is used to determine the flow rate through 
the calibration nozzle. The corresponding voltage from each of the wires 
for each Rotometer setting is then noted. Curves of E versus u can then be 
drawn for each wire. Fig. 3. The derivative 3fcn E/3fcn u Is needed to 
obtain the sensitivity coefficient, that is 

a , MOJ. = J* Ji 

m 3£n u E 3u 

A certain number of points are then taken from each of the two curves 

to obtain A and A . For each velocity a certain D exists, knowing this 
mj m2 m 

a graph can be drawn of D m against velocity and can be approximated by a 
straight line (Fig. 6). From this, any D m is known within the ranqe of 
velocities encountered. 

To obtain the sensitivity to angulation, A y , the prooe is located 
on the x-z plane of the jet (see Fig. 4). The probe is then rotated be- 
tween +10° and -10° incidence angle (#) to the free stream flow direction 
in increments of 2°. This Is repeated at the same velocities at which 



correlator output 

B T r V 
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the D m 's are calculated. A sample of these readings are shown ir Fig. b. 

It is noticed that these curves do not intersect at the zero degree point 
ss one might expect, due to the Inexact matching of the wires - slight 
differences In Individual wire output for Identical flow conditions. From 
the plots of E against ♦, the derivative 3£n E/3$ Is obtained, hence the 
sensitivity coefficient Is also obtained. 

A . 3*Ll. I 

v h l a* 

From each of these graphs A and A u can be calculated, therefore 

V i V 2 

leading to B^. This again Is plotted against velocity and can be approxi- 
mated by a straight line as shown in Fig. 6. 

3.3 F re e Jet Measurements - Proof of Principl e 

Preliminary measurements were made in a free jet for which much pub- 
lished data exist. This Is necessary to validate the experimental pro- 
cedure and data reduction discussed earlier, as little published data con- 
cerning confined jet fluctuating flow data using the present method are 
available. As mentioned previously, measurements were made ^f mean velocity 
and turbulence intensities in the self preserving region of the jet i.e. 
at an axial distance, x/d, of 20. The-e results were then compared with 
those of Con s in. 7 Corrsin 7 used a one inch diameter jet issuing into 
still air at a velocity of 10 m/s while the present jet consisted of a 
contoured nozzle with an area contraction ratio of 16:1 and an exit velo- 
city of 14 m/s. 

Initially, traverses for mean velocity were made across the jet at 
x/d equal to 20 to indicate if the jet used in this study is similar to 
the one used by Corrsin. 7 Figure 7 *ndica f es that the present jet is very 
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compa table with the previous researcher's Jet. Ihe differences that occur 
at the large radial distances can be expected due to the diffi . lties in 
measuring very low velocities. 

Presented in Tigs. 8 and 9 are measurements of u' /u and v* /u 

rms m rms m 

at x/d * 20. As can he seen all the jet measurements have been non-dimen- 
sional 1/ed with regard to the station maximum axial velocity, i.e. on the 
jet centerline. This is because non-dimensional izing with the jet exit 
velocity would be very insignit icant at such a distance downstream from the 
nozzle exit. As can be seen from the plots fair agreement with Corrsin 7 
was achieved, expecially for u ' c /u. The v' _/u m values measured here 
tend to be slightly lower than those presented hy Corrsin 7 but this again 
can be explained by the difficulty in measuring turbulent, low speed flows. 

3.4 Confined Jet Measurements 

before any measurements wore made in the test facility the cros ed 
hot-wire probes had to be calibrated in the free jet facility as described 
in section 2.2. This was tc study the angular sensitivity of the probe 
and to ensure that no 'drifting' of the electronics had occurred. The 
resistances of the hot-wires were also checked after the final measurements 
for the same reason. 

The measurements of mean and fluctuating components of velocity were 
made in an isothermal model of a can combustor The expansion ratio D/d 
was equal to 2 with an expansion angle, a, of 90°. The measurements were 
made at axial locations in the test section of x/D of 0, 0.5, 1, 1.5, 2 
and 2.5. The fluctuating components of the flow were normalized with the 
jet nozzle exit velocity, u o> This quantity was measured prior to the 
commencement of any test. 
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Measurements of u* and v' were achieved by placing the crossed 

rms rms 

wires of the probe in the x-y plane (see Fig. 4) at the required x/D loca- 
tion and Instantaneously adding and subtracting the rms voltages from 
each of the wires. The component of the fluctuating velocity was 

measured by rotating the crossed hot-wires through 90° until they were on 
the x-z plane. The subtraction of the two instantaneous nms voltages is 

then used in the calculation of w' ,,/u . 

rms o 

The probe was traversed radially across the test facility by means 
of a probe drive built at OSU. This probe drive made It possible to 
make measurements at 18 positions tadially across the test section. 
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4. CONHNED JET RISULTS 


The results presented here are measured values of the mean axial velo- 
city, axial, radial, and tangential velocity fluctuations found In isothermal 
airflows of axisyrnnetric combustor geometries. As noted earl ler the diameter 
expansion ration D/d * 2, inlet Reynolds number Re rf ■ 1.26 x 10 s and wall ex- 
pansion angle, a, is 90°. Measurements were taken at axial locations of x/D 
equal to 0, 0.5, 1.0, 2.0 and 2.5. 

Figures 10 (a) - (f) show measured values of the time-mean axial velo- 
city at all of the axial locations. The axial location of x/D of zero is 
actually a distance of approximately 5 mm between the wires and the enlarge- 
ment face to insure that the wires were not damaged. The measurements at 
this location show that the shear layer is very thin, in the region of 3 mm. 
Comparisons have been made with Chaturvedi 10 at all axial locations except 
at positions of x/D equal to 0 and 2.5. Chaturvedi 10 did not take measure- 
ments at these locations. Both sets of axial velocities show good agreement 
except in the region of the recirculation zone. This is due to the hot-wire 
not being unable to sense flow direction, but an increase in flow velocity was 
detected as the probe was traversed further across the recirculation bubble. 
Therefore, Indicating that a flow r°versal did take place. 

The plots depicting u' rms /u 0 are shown in Fig. 11 (a) - (f). As can be 

seen, good comparisons can be noted for most of the measured values with 

those of Chaturvedi. The curves of radial velocity fluctuations are shown 

in Fig. 12 (a) - (f). The two studies are in fair agreement except in 

regions of high turbu’ence where the present study gives higher values 

cotrpared to Chaturvedi. 10 Fig. 12 (a) (v' ,/u at x/D ■ 0) indicates 

rms o 

that the radial velocity fluctuations are damped due to the presence of 
the enlargement wall. 
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The final curves. Fig. 13 (a) - (f) show the tangential velocity 

fluctuations (w* rms / u 0 ) • Chaturvedi 1 °only made preliminary measurements 

of these fluctuations and concluded that v' * w' . The figures shown here 

tend to support this fact with very similar peak and centerline turbulence 

Intensities between v' ,/u^ ana w ' Vu . Fig. 13(a) again shows the 

rms o rmsT o 

damping effect of the wall, which lends to provide further proof of 
Chaturvedi 's l 0 observation. The snail discrepencies that do occur between 
the present results and those of Chaturvedi 1 0 could be caused by slight 
drifting of the electronics and differing nozzle geometries. 

The results shown here show the high level of turbulence in zones of 
separation, the conditions of nonhomogenl ty and the variation of turbu- 
lence characteristics as the flow proceeds. Chaturvedi* 0 noted that center- 
line values of turbulence intensities reached a maximum at an axial loca- 
tion of approximately x/D * 2.5. This result agrees favorably with that or 
the present study. Also shown on Figs. 10-12 are single-wire six-orientation 
measurements , recently obtained in the same test facility at Oklahoma State 
University .* k Again the agreement is generally satisfactory. 
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5. CO NCLUS IONS 

Application of the crossed hot-wire anemometry technique to turbulent 
nonswirling flows has been discussed, together with the calibration and 
data reduction procedures. Great care has to be taken during calibration 
of the probe and, in conjunction with this, 'drifting' of the electronics 
has to be kept to a minimun. The method has oeen applied initially to free 
jet flows in order to justify the technique; results are in good agreement 
with Corrsln. 7 The main application of the method is to axisyrmetric confined 
jets in the Oklahoma State University test facility. Results are given for 
the nonswirling flow in t.he sudden expansion configuration. Measurements of 
time-mean and fluctuating velocities are presented at several axial stations 
and compared with those of Chaturvedi 10 and Janjua." Good agreement is to 
be noted. 
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7 . NOME NCLATURE 

A^ axial velocity fluctuation sensitivity coefficient 
A y radial velocity fluctuation sensitivity coefficient 

»♦ < E ' A v, * E ' V 

D chamber diameter 

d nozzle exit diameter 

< E > V * E > V 

E mean voltage 

e' voltage fluctuation 

r radial distance 

r 0 free jet nozzle radius 

R, hot-wire resistance 

w 

u time-mean axial velocity 

u„, station maximum velocity 

m J 

u„ nozzle exit velocity 

o 

u' „ rms axial velocity fluctuation 

rms 

v' „ rms radial velocity fluctuation 

rms 

w 'rms 0115 tan 9 er, tial velocity fluctuation 

lT' v ' rx component of Reynolds stress tensor 

x axial distance 

x,r,0 axial, radial, polar co-ordinates 

x,y,z Cartesian coordinates 

a wall expansion angle 

<t> yaw angle of cross wire probe relative to jet centerline 
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8. FIGURE LIST 


Fig. 1. Schematic of confined Jet facility 

Fig. 2. Schematic of Instrumentation layout 

Fig. 3. Calibration plot of E versus u 

Fig. 4. Coordinate system of free and confined jet 

Fig. 5. Calibration plot for sensitivity to angulation 

Fig. 6. Parameters D and B n versus u 

mo 

Fig. 7. Radial distribution of time-mean axlai velocity at x/d * 20 In 
nonswirling free jet [ Ref. 7] 

Fig. 8. Radial distribution of axial velocity fluctuation at x/d * 20 in 
nonswirling free jet [ Ref. 7] 

Fig. 9. Radial distribution of radial velocity fluctuation at x/d * 20 in 
nonswirling free jet [ Ref. 7] 

Fig. 10. Time-mean axial velocity profiles In ncnswirling confined jet 
[a — o present study, Ref. 10, Ref. 11] 

Fig. 11. Axial velocity fluctuation profiles in nonswirling confined jet 
[o o present study, Ref. 10, Ref. 11] 

Fig. 12. Radial velocity fluctuation profiles in nonswirling confined jet 
[o — o present study, Ref. 10, Ref. 11] 

Fig. 13. Tangentia 1 velocity fluctuation profiles in nonswirling confined 
jet 
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Fig. 2 Schematic of instrumentation layout 
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Fig. 5 Calibration plot for sensitivity for angulation 
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Radial distribution of time-mean axial velocitj 
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Fig. 8 Radial distribution of axial velocity 

fluctuation at x/d = 20. [ Ref. 7] 
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A SIMPl F FlfilTF DDMHfNlF PkOMMIRt 
FOR THE VORTfX OOKTROILFD 

A. A. Rusnalna* and D. li. Lllley** 
Oklahissa State University, Stillwater, Ok la. 


Abat ran 

The effect of hired on ahort vortex controlled 
dtffusera la considered numerically. the tr.inatanl 
Nav let -Stoke* equation* arr solved via their asao- 
clrlril llnlta difference aqua! Iona, dlrartly In 
taima .»l tha primitive pressure-velocity variables, 
•ind I ha technique la baaad on (ha two-dimensional 
Los Alamoa SHIA prediction proiadura. A ahort alm- 
pla coaputar coda haa baan davalopad, using a l.ual- 
nar Mow simulation with 'fra* allp' or 'no allp' 
wall boundary condition*. rha alapllltad coda re - 
piaaanta a baatc tool to which user-oriented roo- 
plaxltlaa and aophlat Icat Iona can aaallv ha addad 
a* rrqullad. 


1. Introduction 
t.l fha Phenomenon Con a Ida rad 

Recant emphasis ha* baan placed on tha require- 
ment tor a ahort efficient dllluaer to be placed 
between tha compreaaor and comhuator component a of 
the aircraft gas turbine engine. The vortex- 
controlled dllfuaet VCD employ* bleed at (he atep 
location of a dump diffuser, ao aa to obtain a low- 
preaaure loaa In a ahort length. Figure l ahowa a 
achematlc of the VCD In it* aimpleat form. Here, 
the view may be conaldered aa a allce ot an annular 
path with Inner and outer bleed* or aa an axlavm- 
tretrlc situation with cl rcuml erent lal bleed. In 
practice, the I low diffuse* 1 tom the primary duct 
Into an abruptly-enlarged aacondarv duct. Ini r - 
actlon between the VCD bleed and tha main I low la 
atabll/ed via a vortex fence, alao ahown In the 
figure. Without bleed, there la )et flow with vary 
little lateral dllfualon - a long duct la nerd to 
obtain uniform air-flow propertle* at the dlffuaer 
exit. Conversely, with sufficient bleed, the main 
flow diffuse* rapidly and Mila tha secondary duct - 
the corner recirculation zone 1* ccnsidrrabl v 
shorter. 



Fig. 1 Schematic of tha VCD teimrtrv 1 

• Graduate Student, School of Mechanical and A. /o- 
spaca Fnglneettig, Member AIAA. 

•* Associate Professor, School of Mechanical and 
Aerospace Engineer Ins, Aaaoefate Fellow AIAA. 

t .>e«M|Xi \at.M ■« ImiilnU of lensisisr 
AwroesaMc*. la* ISR1 All n**i» >» »*»»d 


1.2 Previous St udie s 

Farlv Ideas of the VCD placed cusps In tha dlf 
fuser wall to locate vortices rotating preferential 
ly In tha direction of flow, ao aa to reduce the 
boundary layer shear rtrasa experienced by tha flow 
In (ha region of high adverse pressure gradients. 1 
(Inly limited success waa obtained because ol diffi- 
culty In retaining a stable vortex ayatae. Recent 
extensive studies at Cranfteld Institute ol Techno- 
logy overcame this problem by featuring hired off 
directly lr<*a the vortex. *•*•' It waa alao found 
that (he quant It of bleed required could be kept 
at teaaonable levels by locating a fence dlrettly 
downstream of the vortex, thereby creating a par- 
tially enclosed vortex chamber. Figure 2 shows the 
mechanism deduced from exper tmental data, where a 
layer of high turbulence forma between the drawn 
off stream tube 's' and the diffusing stream tube 
'b'. The mainstream flow la further aided on Its 
path by the formation of a found* bubble immediate- 
ly downstream of the lateral fence. 

This work suggested a potential performance 
payoff In practical rnglne appl Icat ions. Detroit 
Diesel Allison Division of General Motors applied 
the above laboratory concepts to a realistic gas 
turbine annular rombustlon system flow path. The 
performance was confirmed to be a strong function 
of several variables Including: 

1, Secondary duct length L/H 

2. Meed quantity R (flow fraction) 

1, VCD area ratio Ag 

4. VCD Inlet flow distortion a 

1. Bleed slut axial and radial gaps 

6. Bleed cavity geometry. 

Other workers considered applications to variable 
cc.snetry combustors*”’ and Inlet flow distribution 

. licit*.* 



1.1 The Computational Problem 

Mathemat leal ly the governing equations are 
ell'ptlc in character and a relaxation method of 
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solution Is appropriate (II solving directly for 
the steady-state) ; they are parabolic (In time) and 
a earthing aethod o( solution la appropriate (II 
solving lor the steady-state via the flee evolution 
of the transient flow process). Consideration la 
given here to a primitive, pressure-velocity varl- 
tlbe, finite difference code which has been develop- 
ed to predict 2-D aria vame trie transient flows. 

The technique uaed here tncorporatee the following) 

(I) A finite difference procedure la uaed In 
which the dependent variables are the 
velocity coaponenta and pressure, formu- 
lated In aalswetrlc cylindrical polar 
coordinates. 

(11) The pressure Is deduced froa the contin- 
uity equation and the latest velocity 
field, using a guess end-correct proce- 
dure for the latter. 

(Ill) The procedure Incorporates displaced grids 
for the two velocity components, which 
are placed between the nodes where pres- 
sure Is stored. 

The code Is designed lor pernons with little 
or no experience in numerical fluid dynamics with 
the purpose of demonstrating that many useful and 
difficult problems can be solved without resorting 
to large, duplicated computer programs. As speci- 
fic complexities can be added easily, the program 
provides a basis for developing many new numerical 
capabilities. The present work has grown out of 
previous Ideas about solving 2-1) flows using so- 
called Los Alamos 1 ®' 1 * and Imperial College 1 *' 1 * 
methods, solving for the primitive pressure-velo- 
city variables, the former treating the transient 
case and the latter treating the steady state case 
directly. 

The present study illustrates the use of a 
short simple computer code, using a laminar flow 
simulation with 'free slip' or 'no slip* wall bound- 
ary conditions. The simplified rode provides an 
alternative to a more sophisticated rode applied to 
a similar dump flowfield with bleed which is dis- 
cussed as one of the gas turbine combustor design 
challenges of the llBO's. 1 * 

1.4 Outli ne of the Con trib ution 

The computat ional procedure described In 
Section 2 is baaed on Los Alamos SOIA 2-D ideas. 1 * 
The technique is for incompressible fluid flow, and 
complexities are restrained to a minimum so as to 
facilitate the use of the rode. User-oriented com- 
plexities can be added as required and details of 
many such extensions may be obtained elsewhere: 
control volunie and flux approach to obtaining the 
FDEs, non-uniform grid system, seml-lmplirit method 
for each time-step, etc. 't 14 * 1 * Predictions for 
the VCD are given in Section 3 showing the influ- 
ence of several design parameters. The final con- 
clusions illustrate that a useful and valuable tool, 
based on conceptually simple ideas. Is now becoming 
available, and on to which any user-oriented com- 
plexities can easily be added. 

2 . SI m u 1 a t Ion and So lu 1 1 on Pr ocedu re 
2.1 The Technljjue 

The computat tonal code solves directly for the 
primitive pressure and velocity variables, and is 
based on the Marker and Cell MAC Los Alamos techni- 


que. This la one of the moat well-known methods to 
solve t 1 me -dependent Incompressible fluid flow pro- 
blems; Its conceptual simplicity Is one of the main 
attributes. A 2-D version Is described In Aef . 12, 
upon which the present work la baaed. In the pre- 
sent work, and Eulerlan finite difference formula- 
tion la used with pressure and velocity as the main 
dependent variables. In addition, the velocity 
components are positioned between the nodes where 
pressure and other variables are stored. At each 
step, the time-advanced expressions far u and v are 
substituted Into the finite difference form of the 
continuity equation for sach cell, and the guess 
and correct Iterative process on pressure and velo- 
city corrections Is dona until the continuity equa- 
tions Is sufficiently well satisfied. This section 
continues with a full dlscuaalon of the partial 
differential equations, mesh system, finite differ- 
ence equations, stability and accuracy, and boun- 
dary condition Insertion. 

2_,2 The P artia l Diff e renti al Equat ions 

The governing partial differential equations 
PDFs of the Incompressible flowfield are the con- 
tinuity equation, and conservation of momentum In 
the x (• radial) and y (• axial) directions. These 
may be taken in axisvanetric cylindrical coordinates 
in conservative form 1 *' 1 ® as: 



The above simulation lndl rates that the pro- 
blem Is simply treated as a laminar viscous problem 
with unrealistic 'Tree slip' or 'no Blip' wall 
boundary conditions. In this manner of formulation 
much of the flow Is controlled by the Inlet and 
outlet condTtions, together with the confining con- 
straint. As recommended In Kef. 12, when computa- 
tions are made in a large chamber, it Is often a 
complicated and lengthy task to resolve thin boun- 
dary ivers along the confining walls. But in many 
cases, the presence of these Is unimportant and 
free slip boundary conditions can be justified as a 
good approximation. On the question of turbulence 
modeling, it is known that the time-mean behavior 
of a certain turbulent flows may sometlsies be simu- 
lated via a constant large viscosit y.*®'* 1 This Is 
the approach taken here and predictions have been 
made with a range of values given to the constant 
viscosity V with minor effect on the flowfield pro- 
duced. In view of this minor effect of the visco- 
sity magnitude, computations given later are obtain- 
ed merely with the laminar viscosity being used. 

The results show that, even with this limitation, 
a useful simulation of the actual problem is ob- 
tained . 
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Th* present simulation and solution plot# Jura 
la a abort simplified approach (approximately )%0 
fortran atafeswntsl with pood capability. Of course, 
othar aura lengthy | approa laata ly 2000 fortran 
atataaantal and aophlat Icatad | nonuni fora grid, la* 
pllclty aolutlon procadura, two-aquatlon turbulence 
aoda l and wall functional approaches ara aval labia. 

2 . ) Tha Ha *h System 


Tha a-y domain la dlvldad Into rectangular 
call dlvtalona, with uniform As and Ay sparlngs. 

In tha caaa of an annular dlffuaar, tha laft hand 
boundaty rapraaanta a symmetry plana with radlua 
»[ • ( r | ♦ r 0 )/2 whara t\ and r„ ara Innar and 
outar primary duct annulua radii. Th 1 a aaaumad 
avmmatry Impltaa that only half of tha pioblam la 
being aolvad. In tho caaa of an axlsyssaet rlc dif- 
fuser, tha laft hand boundary la tha symswtry axla 
with radlua • 0, 

Thla domain la alao rumplcmantad by a layar of 
alia on all aldaa, ao aa to allow aaay almulatlon 
of tha required boundary condltlona (BCs). Thaaa 
flctltloua calla Increase tha total numbar of calla 
In each direction, figure 1 Uluatratea the total 
aaah arrangement , while fig. 4 portray# a alngle 
call and ahuwa tha location of each field varlablea 
p, u and v relative to thla (l,J)-cell. Tha pres- 
aura p la located at tha canter of each call, while 
tha radial and axial valocttlea (u In tha radial 
x-dlrectlon and v !n tha axial y-dlrectlon) ara on 
tha right and top boundarlea, respectively. Thua 
normal velocities lie directly on the phyalcal 
boundarlea of tha aolutlon domain, while the tangen- 
tial veiocltlea and pressure arc dttplaced half a 
cell interval Inalde the flowfleld. In thla wav 
the exterior flctltloua calla are particularly i un- 
vent ent whan applying the boundary condition# BCa. 

.'.4 The Fini te J)l_f f ereiu «• I <iujt 1 ons 

In formulating finite diffarence repreaenta- 
tlona of tha governing partial differential equa- 
tlona PDEs, the uaual intuitive catlmatea of ore- 
aided flrat -derivative#, cantered f lrst-derlvat ' ves 
and centered second- derive! Ivaa are uaed. Super- 
acrlpta n and (blank) are uaed to denote valuea at 
time- level t and t ♦ At , respectively. Portrayed 
now are equatlona enabling one auch forward t Ime- 
atep to be accomplinhed. Thua, atartlng (torn Ini- 
tial field valuea throughout the domain of Intereat, 
a time-march proceau la used ao aa to advance to- 
ward the final ateadv-atate aolutlon, which la 
uaually of apcclal Intereat aa oppoaed to tha en 
route calculat Inna. 

In Eq. (1) the time derivatives are approxi- 
mated by forward one-alJed derlvatlvea; av«t apa- 
tlal derlvatlvea are approximat'd by cential differ- 
ences baaed on valuea at Claw-level t. Special 
lechnlques are required In computational fluid dyna- 
mics, however. In the representation of the convec- 
tion terma, and a certain amount of upstream dif- 
ferencing la required. The difference equations 
representing the conservation of momentum Eq. (1) 

In directions x and y may be written and used to 
accomplish one forward time-step as follows: 



Fig. 1 General grid system covering the 
flow domain. 
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Fig. 4 Arrangement of finite difference 
variables In a typical cell. 
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vtirir the * l«nt on ilit right of itch of these 
tqutl lont tit defined in Appendix A. Th* Coot f 1- 
rient a in the** equations Appendix A) io • 

amtltnl taking a value be teen 0 and 1, and ao giv- 
ing l ha daaltad a noun I ol upatraaa (donor rail) 
dlffarancing fn I ha ronvaclion terms. 11 A valua of 
0 glvas merely rantral diffaranclng as in tht ori- 
ginal MAC rodaa 1 * > * 1 and nunarlcal lnaiability pro- 
blana ariaai a valua of 1 givaa tha full upatraaa 
or donor rail fora which, though Introducing rnoit, 
la atabla prnvldad tha fluid ia not allowed to paaa 
through nor * than ona rail in ona tlaa-atap, aaa 
tha atablllty rritarion of tq, (b) latar. 

Although Eq. (2) arronpliahaa ona forward tlaa- 
atap baaad on ronaarvatlon of momentum prlnriplaa, 
tha nawly ralrulatrd valor It la* will not, in gen- 
aral, aatlafy tha continuity raqulraaant, aa ex- 
pressed by tha rantral I' ilia difference fora of 
tha rontinulty aquation (aaa tha firat of Eq, (l).t 


A* <U 1 J * u l-!.) > 4 Sy <V 1J ’ V I,J-1 ) 


x 

2x 


(I) 


11 
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ralaaatlon paraaatar w batwaan 1 and 2, a valua 
nerr l.l oftan being optlaal, in ordar to aoaad Uf 
tha cmvetgenc* of tha pnaaura Itaratlon procaaa. 
Convergence to |ha alaady-atal* aolutlon la eatab- 
llahad by taking aanv forward tlaa-atrpa. 


2.) St ability and Accuracy 

Tha ataady atata aolutlon ia avolvad Iron 
taking many forward tlaa-atapa, a rarnaaandat ion 
ng at laaat twira tha avaraga raaldanra tlaa of 
a typlral partlrla paaalng through tba flow domain, 
baaad on a a lap la 1 -D flow analyaia. Arrurary la 
anhanrad by ualng aaall apatlal and tlaa Intarvala, 
at tha aapanaa of larga coaputar tlaa. Whan total 
tlaa raqulraaanta naraaaitala a larga grid alt*, it 
la not poaatblf to raaolva thin boundary layara 
along ronfining walla and fraa-allp BCa for tangen- 
tlal valorlllaa ara aura appropriate than tha no re 
Intuitively rorrecl no-altp BCa. 


Once tha apatlal aubdlviaion ia rhoaan, tha 
tlaa Inrreavnt auat be raatrirtad in two way*. 
Firefly, it must ba laaa than (typlrally equal to 
0.2) to 0.1) t laaa) tha alnlania rail transit 1 1 aw 
taken over all relle: 


At < aln 



<6) 


Secondly, with a nonrrro klnenatir vlscoalty, momrn- 
tua auat not dlffuke aora than approalnataly ona 
rail in ona tlaa-atap, for whirh an *• aata ia: 


V At < 


1 

2 



( 7 ) 


where 


» • « ( ♦ Ah (1-1. 5) 

Tarn* hare are evaluated at tlna-laval t ♦ At . 

Thla inrumpreaalolllty condition la laipoaed by 
Iteratively adluatlng tha rail pressure. That la, 
if the divergence D of a cell ia poaltive (the left 
hand aide of Fq. (1) ia poaltive) there la a net 
nave outflow Iron 'hat rail. Thla la correrted by 
reducing the cell preaaure. If the divergence la 
negative, an increase in rell preaaure la appro- 
priate. 


When the t law-step la ao restricted, the re- 
quired amount of upatream (donor rell) dlffetenc ng 
must hr achieved by chunking a kllghtly larger than 
(typically 1.2 to 1.) time*) the largest of the 
right hand aide member of: 



(B) 


where the maximum la taken over all cell*. If a Is 
chosen to be too large, atablllty ia being achieved 
at the expense of the intorduction of an unneces- 
sarily large amount of diffusion-like truncation 
error* (called numerical aauiothlng). 


When a call preaaurr changes from p to p a Ap, 
the velocity components on th* 4 fares of that rell 
change, given from a linear analysis from Eq . (2), 
by an amount: 


u . • u. . ♦ At Ap/Ax 
1,3 1,3 

Vi.) • “i-i.j * At Ap/A * 

V 1,J ■ v l,j 4 At Ap/Ay 

v i.J-l " v l,j-l - At A,,/Ay 


(4) 


2 . 6 Bounda ry C ondlt ions 

Finite difference equations FEDs simulating 
the PDF.a are aet up and solved by way of a time- 
march proceaa applied to cells within the flow do- 
main of interest. Cells touching the boundary thus 
util ire the value on or just hevond the boundary. 

Interior normal velocity calculation* take the 
aero normal wall values, the given normal inlet or 
outlet value*, or the vet-to-be-determined outlet 
values aa appropriate BCa during their calculations 


Substitution of these in Eq. (1) yields the 
amount of correction to p required as: 

Ap - - D/(2 At (1/Ax* ♦ 1/Ay*)) (1) 

where D is the current (nonzero) /alue of the left 
hand aide of Eq. (3). 

At each time-step, pressure and velocity up 
date iteration continues until the Da of all the 
cells are less than some prescribed small poaltive 
quantity C. Th* Eq. (4) are applied with an over- 


interior tanqentiaJ velocity calculations use 
‘.he fictitious values which are placed in the sur- 
rounding layer of coaplementary cells. Specifica- 
tion of these la after each tiaw-step and after 
each sweep of the cel 'a during the pressure itera- 
tion. With a coarse gild, frm-slip BCa are appro- 
priate for tangential velocities, and external 
values are set equal to their associated immediate- 
ly interior values. On the other hand, with a fine 
grid computing through the boundary layer, no- a l ip 
BCa are appropriate for tangential velocities and 


4 


external valura arc equal to the negative of 

(halt aas-c loted Immediately tntarlor valuaa. 

Specif lcat Ion ol normal valocitlaa at an tu t - 
flow boundary of tan poses a problem aa It can hava 
detrimental upatrr.ua influanca. Ona al|)it merely 
Impure t ha zero-normal gradient or cont inujt iva 
condition and aat thaaa valuaa aqual to thalr Im- 
mediately upstream valuaa. 1 * When primary Intaraat 
la ba ln| focuaad on tha final ataady-atata ablution. 
It haa baan found 1 ’ that a aultabla constant ray ba 
addad to each such extrapolated value, with advant- 
■pe to tha rapidity of converRence. This constant 
value la chosen so aa to nakr tha total outlet f lus 
aqual to the total Inlet flus, thua ensuring the 
requl rrnrnt of a aacroacroplc naaa balance. Outlet 
boundary specification la Imposed only after each 
time-step aa computed via Fq. (2) and not after each 
pass through tha me ah during the pressure Iteration. 

On an axis or plane of symmetry the usual tero 
normal velocity and free-sllp axial velocity speci- 
fication are applicable. 


I I’r ed U; 1 1 ona and Di scussion 


1.1 Diffuser Ef (act lveness 

Experimental tests 1 * 1 have Included both non- 
uniform Inlet profiles and varying quantities of 
bleed. Modifications to the usual ef fact lveress 
parameter are needed to characterise the situation 
correctly. Flow distortion la characterized by the 
flow distortion parameter, a. This Is defined as 
the kinetic energy flux of the given profile divided 
by the kinetic energy flux of the same flow quanti- 
ty with a flat profile. Values larger than unity 
Indicate flow distortion. The following equation 
may be used ro compute a: 

« * 


where the Integrations are taken across the flow 
passage. Tills parameter Is of fundamental impor- 
tance since It establishes the ideal diffuser static- 
pressure rise coefficient , which for Incompressible 
flow t s ; 


C - a i - at 
P 1 
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field residence times later, corresponding to 100 
time-steps, and requiring about 70 a of IBM 3)0/168 
CP time. The iM-aber of pressure-velocity Iterations 
required par time-step depenos on tha Inlet nonunl- 
formltlea and bleed magnitudes. Typically, It Is 
Initially about 150 and decreases to 1 as the final 
steady state solution Is approached. 


Figure 5 shows predict lone of the velocity 
vector patterns In the xy-domaln, with OX and 16X 
bleed In Parts (a) and (b), respectively. Notice 
that with no bleed, there te a rather large corner 




Fig. 5 Predicted velocity vector patterns 
with: 

(a) OX bleed 

(b) 16X bleed 

|— recirculation bubble boundary* 


The diffuser effectiveness Is the actual diffuser 
static pressue rise, divided by the ideal static 
pressure rise. Typically, the Inlet flow to the 
diffuser is highly distorted with at - 1.22 and 
Ug/ugy * 1.2 (the maximum axial velocity is near 
the center), while the exit flow Is far more uni- 
form with tfj and u a /u av taken near unity. 

1.2 Predictions 

To demonstrate the capabilities of the solution 
scheme Just described, predictions are now made and 
discussed for the case of the sudden expansion an- 
nular VCD with area ratio Ag - 3, 2 and 1.5. The 
Inflow is taken to be highly distorted with at *1.22 
as in the experimental data used for cosqiat Ison. 1 

Predictions are first discussed for the case 
Ag - 2. The flow domain Is overlaid with 10 x 28 
Internal cell divisions in the x- and y-dlrections 
respectively, with aspect ratio Ay/Ax • 4. Tran- 
sient computations are continued until steady stale 
conditions establish themsu 1 ves , roughly !»■: flow- 


recirculation zone provoked by the sudden expansion. 
It Is some distance downstream before the axial 
velccity profiles becoaw more uniform. On the other 
hand, with 16X bleed, the Inflow Is encouraged to 
expand Into the larger cross-sectional area and only 
a relatively small corner vortex region Is to be 
seen. Axial velocity profiles c ecome more uniform 
much more quickly In the downstream direction than 
was the case with no bleed. The dashec line on 
theae vector velocity plots Is approximately the 
zero streamline, which enclosed the comer recircu- 
lation bubble. 

Vortex controlled diffuser performance predic- 
tion le Illustrated In Figs. 6-8 for the are ex- 
pansion ratios Ag • 3, 2 and 1.5 respectively. In 
each case, VCD effectiveness le pic, ted as a func- 
tion of downstream distance L/H, measured from the 
end of the primary duct. Each of the figures has 
three lines drawn, representing thf cases of OX, 

81 and 12X bleed. The open symbol! and solid lines 
are predicted; the closed aynlx-is .ire corresponding 
expei latent a 1 evidence. The agreerent la very 




Mg. t> VCD pet I muncr vllh Mf> i*t In Ag • 1 
and Inlrl I low dlalortlon a • 1.27. 



DUCT l/H 

Kl||, 7 VCD performance with ttra ratio A. • J 
and Inlrl flow dlalortlon a * 1.7*. 



DUCT L/h 


FI*. « VCD performance with arra tatlo 

K * 1.5 and lnlrt flow dlalortlon 

O - 1.22. 

aatlalartory In virw of thr almpllrltv of thr «■- 
putrr tod* being uard for thr a I anil at Ion and thr 
general trrnda arr cirarly rvldrnt. 

An rarltrr atudy 1 aaaerta |but dura net clearly 


do< uarent | that when rf I act tvrnraa la p lot t rd aga Inal 
thr prrcrntagr ol main flow tha> la bird from thr 
dlffuaet, llirlr la a ‘mtriiiaua Meal nyulirari I'. 
Rlrrd prtirntagra hrlow thla valor arr Ineffnrtvv 
and, aa thla Irvrl la rrarhrd, thnr la a ratlin 
auddrn Inctraar In prtfotmanir. Abovr thla valor 
only a alight forthrr Inptovrnrnt la obtained Thr 
prrarnl pirdlrtrd rraolta do not ahow thla abrupt 
trrnd, and thr onarl of prrforBanrr Inctraar la lat 
lraa dlaaatlr. Thrrr la a ront Inooua lapiovrarnt 
In prrlonaancr aa thr blrrd pneantagr la Increase-*. 
Figure *i ahowa thr drtalle for thr area ratio r<|oal 
to twti caar, with corvrd of prrfotaancr vrtaua 
blrrd. Thr two llnra rrptrarnt pnforautner bring 
rharai trr lird by valura at local Icna l./H • 7 and I. 



Ktg. ** VCD performance with >ica ratio Ag • 2 
and Inlrl flow I'.latortlon n • 1.22. 

A. Cone lua Iona 

A atmplr prediction procedure for auddrn r» 
panalon Im ompreaatblr I Ioum haa been developed end 
applied to thr vortex controlled illlluart. H.raeil 
on thr 1 ia Alamo* 2-D SHI A Idraa, thr ttanair>nt 
Navler-Stoke* equation* ol an Ineomptraxlblr fluid 
air aolvrd via their aaaoclatrd finite dlffrrrncr 
equal Iona direct Iv In trraa of primitive prraauic- 
veloclty variable*. Thr technique ia atmpltftrd to 
ladlltatr lta uar bv prrauna with little or no ex- 
perience In comput at tonal fluid dvnamira. and boun- 
dary rondltlona, tor example, air particularly aim 
plr to aprclfy. Thr code thua rrprrarnta a haait 
tool, to which uart -or tented complrxltlra and ao- 
phrat teat lo'ia can easily be added aa required. Thr 
rraulta pcaented and dlacuaard arr extremely uae* 
ful to thr dralgner of vortex confrolltd dlfluarra 
for pte-combuator appl leal Iona . Predicted rraulta 
confirm that ef feet tvrnraa Incrraara with Increaara 
In dart lrngth end bleed f low rate. 

Other development* and application* of thla 
finite difference technique currently In progim* 
at Oklahoma State Dnlveraltv Include: Inclualon 

of awlrl and thr two-equation k-i turbulent.' model 
Into thr preaent 2-D code, and laminar and turbulent 
3-D code* for application to local dratrat If leal Ion 
and dilution air croaa-|rta Into main attrama. 
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Append I x A 

The finite dllfe.ence form of the Navler- 

Stokes equal Ions *s given as Fq. (7) with donor cell 

differencing In sixnc of the terms. For complete- 
ness, the appropriate expressions are appended. 
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Nomencla ture 

A Flow atea. »' 

A h Dtftuaer arra expanalon rrtio 

R Dill user Mrrd (low tia t Ion 

C pi ldral dlffuarr atatlc preaaure riar 

corf flclrnt 

1) Plvrritrnre o( velocity field, a" 1 

H Channel height (width), m 

t. J Meah point 

L Channel length, B 

p Ratio o' prraaure to ronatant denaity, 

Na/kg 

t Tlar coordinate, a 

u, v Velocity componenta in x( r),v 

coordinate dtrectlona, a/a 
x,y Cylindrical polar coordinatea, a 

Ap Small finite p Increment, Na/kg 

At Small finite t lae lncrraent, a 

Ax.Av Small finite radial and axial dlatancra.a 

ci Donor cell upatrram differencing para- 

meter, ' ’ ow dlatortlon factor 
C S aa 1 . aitlve nuaber 

V Ktnra .tic vlaroalty, a : /a 

p denaity, kg/a’ 

u Over-relaxation parameter 

V I Laplaclan operator 


Subacrlpta 

av Station average 

rh Chamber 

1,1 Meah point Indexea 

i Inner 

a Station aa-daua 

o Outer 


♦ At 
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A BASIC CODE FOR THE PREDICTION OF TRANSIENT 
THREE-DIMENSIONAL TURBULENT FLOWFIELDS 
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ABSTRACT 

A primitive pressure-velocity finite difference code has been 
developed to predict transient three-dimensional turbulent flow. 

The code is a simplified yet effective prediction procedure for use 
by persons with little experience in computational fluid dynamics, 
and into which user-oriented complexities can easily be added. The 
method is based on the transient two-dimensional Los Alamos SOLA 
prediction technqlue for laminar flows. Turbulence is simulated 
by means the two-equation k-e turbulence model; species diffusion 
and buoyancy are also lnnuaea. T wo applications of the code are 
presented to local destratification near the release structure of 
a reservoir and to the deflection of a jet entering normally into 
a uniform cross-flow. Predicted results exhibit good agreement 
with experimental data, showing that a useful characterization of 
fully three-dimensional flows is now available. 


Paper for presentation at the Symposium on 3-D Turbulent Shear Flows 
to be held at the ASME Fluids Engineering Division 1982 Spring Meeting 
in St. Louis, Missouri, June 7-10, 1982. 
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k 

P 

Kinetic energy at point P 

C 1 

k-C model i'onalant 

K 

Conatant In the log-law 

C 2 

k-r model constant 

K« 

Nondlaenslonal propeller distance, nornallt 

S 

k-r model constant 


ed with respect to H. 

C M 

k-r model conn (ant 

L 

Propeller depth below the surface, a 

D 

P 

Propeller diameter, ■ 

( 

Length scale of turbulence, a 

D* 

Nondlaenslonal propeller diameter, C*“D^/H 

l* 

Nondlaenslonal propeller depth, L* ■ L/H 

or 

Dilution factor, DF - (P f -Pj)/(P.-Pj) 

"l 

Hass fraction of fluid of density pj 

l 

Constant in the log- law 


(epl llanlon) 

Frd 

Demtmet rlc Froude number. 

"2 

Hass fraction of fluid of density Pj 


Frd - v (| /|,(Ap/p 1 )H| l/J 


(hvpol Imnton) 

* 

Gravlt tonal acceleration, a/i 1 

P 

Tlae-aean static pre»sure deviation fraa 

C 

Generation rate of turbulence energy. 




J/kg-s 

P 

Typical point in flowfteld 

H 

Total depth of model, ■ 

9 P 

Propeller flow rate, kg/a 

k 

Kinetic energy of turbulence per unit aass. 

9 r 

Release flow rate, kg/a 


J/kg 

Q* 

Nondlaenslonal flow rate, Q* " Q /Q 

P f 



R 

Jet-to-crossf low velocity ratio, R - Uj/u^ 


• 

r 

o 

Deflected Jet radius, a 

• 



Source tera for k 



S c 

Source tera for c 


• 

t 

Tlae , a 
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u.v.w 


V 

9 

V 
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l.T.I 


V 


c 

u 

V 
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T lM-win velocity coaponrnd In 1,1,1 
coordinate directions, a/s 
Propeller velocity , ■/• 

Coaponen ( of velocity parallel to the wall 
In tha near wall region. m/m 
Cartealan coordinates, ■ 

Nornal distance froai nearby wall to point 

r. m 

Thermocllne local , jn (measured (roe tha 
bottom?, ■ 

Mondlnenslonal thermocllne location, 

V - v 

Turbulent energy dissipation rate, m’/s* 

Absolute viscosity. Pa. a 

Turbulent viscosity, V,.m 

Klnenatlc viscosity, m l /t 

Eplllmnlon or top fluid density, kg/a* 

Hypol 1 union or bottoa fluid density, kg/a* 

Release fluid density, kg/a* 

Schmidt number for k 
Schmidt nuaber for C 
Schmidt nuaber 

Turbulent stress tensor, N/a 1 
General dependent variable 


INTRODUCTION 


Existing 3-D turbulent flowflelds prediction 
codes, although Incorporating many complexities and 
being efficient in their solution algorithms, pre- 
sent a aajor struggle to the code user who la faced 
with the task of understanding, amending, and utilis- 
ing the available codes (1). There la a need for 
simplified techniques for persons with little or no 
experience in computational fluid dynamics CFD, and 
Into which user-oriented complexities can easily be 
added. This paper presents a method which. In a 
computer program of approximately 430 statements, 
solves the fully 3-D time-dependent turbulent flow 
equations in cartesian coordinates. Turbulence la 
slsHilated by means of tl>c two-equation k-C turbulence 
model (2), and species diffusion and buoyancy effects 
are properly simulated. The computational method la 
based on the transient 2-D Los Alamos SOLA prediction 
technique (3) (for laminar flows], which is a finite 
difference scheme based on the marker and cell method. 

Conalderat ton Is given here to a primitive, 
pressure-velocity variable, finite difference code 
which has been developed to predict turbulent 3-D 


transient flows. The technique uaed here Incorporetee 
the following: 

1. A finite difference procedure la used In 
which the dependent variables are the velocity com- 
ponents and pressure, formulated in cartesian coordi- 
nates. 

2. The pressure le deduced from the continuity 
equation and the latent velocity field, using a gueaa- 
and-correct procedure for the latter. 

3. The procedure Incorporates displaced grids 
for the three velocity components, which are placed 
between the nodes where pressure and other variables 
are stored. 

4. The code la designed for persons with little 
or no experience in numerical fluid dynamics with tha 
purpose of demonstrating that many useful and diffi- 
cult problems can be solved without resorting to large, 
complicated computer programs. As specific complexi- 
ties can be added easily, the program provides a baala 
for developing many new numerical capabilities. The 
present code incorporates the two-equation k-C turbu- 
lence model (4). 

The paper presents recent work In the finite dif- 
ference solution, via a primitive pressure-velocity 
code, of 3-D recirculating transient turbulent flow- 
fields. The technique Is for incompressible flow, and 
the complexities are restrained to a minimum so as to 
facilitate the use of the code and user-oriented com- 
plexities can be added as required. Two applications 
of tha code are presented. The first Is prediction of 
local destratification of reservoirs (3,6) which In- 
cludes species diffusion and buoyancy forces, and the 
second application la prediction of the velocity field 
generated by a deflected turbulent Jet, injected nor- 
mally into a uniform cross-flow (7). 

THE PREDICTION PROCEDURE 

Covernlng Equations 

Equations representing conservation of mass, 
momentum, species diffusion, turbulent kinetic energy 
and lta dissipation rate are taken in conservative 
form in Cartesian coordinates (6) as: 


*1 

*r° 

I? 4 £[ (u i^ ■ 8 ? + y 


a) 


where 4 stands for any of the dependent variables u, 
v, v, m , k and C. The first three of these are time- 
mean velocity components In the x. directions 
(x, y and t|, respectively. A buoyancy term occurs In 
the vertical (-direction equation (6). The local maas 
fraction mi of fluid of density Pi allows tha local 
density to be found us'.ng 


•I ♦ mi ■ 1 
p - mi pi ♦ mi p« 


(2) 


Density changes are alight and allow lncompresslbla 
equations to be used, with the addition of an upward 
buoyancy force. The viscosity Is calculated from tha 
standard k-C turbulence model (4) 

M t - Cypk*/c (3) 

where k la the turbulence energy and C [ - Cjk’^ll la 
the energy dissipation rate, both of these being ob- 
tained from solution of partial differential equations. 
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specified a priori 


The aquation* differ primarily In thalr Schmidt nua- 
bara tu and (Inal 'source' terms Sa, aa lndlcatad In 
Tabla I, Tha rata of generation ol k by tha action 
o( velocity gradient a la taken ao 


C 



( 4 ) 


Turbulence and other conatanta appearing are glvan 
tha uaual tecoaaaandad valueat Cq • 1.0, C • 0.09, 
Cj ■ 1.44, Cj - 1.92. Thaaa have bean uirH In a 
wide variety of turbulent flow altuatlone and exhi- 
blt good pradlctlva capability (4). 


TABLE 1 


SCHMIDT NUMB IRS AND SOURCE TERMS 
IN THE GENERAL EQUATION 
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The Crld Syaten 

The rectangular 3-D region to be conaldered la 
divided into rectangular cell dlvialona, with ’<ni- 
form Ax, Ay and Az spaclngs. Thla solution donaln la 
compleaented by a layer of cells on all aldea, ao aa 
to allow easy alaulation of the required boundary 
conditions. These fictitious cells increase tha 
total nuober of cells in each direction. Figure 1 
illustrates the total nesh arrangement, while Figure 
2 portrays a single cell and shows the location of 
each field variable p, u, v, w and relative to 
this (I,J,K)-cell. Notice that pressure p and 
are cell centered, as also are the turbulence quanti- 
ties k and £. Velocities are located on the faces 
of the cell such that appropriate pressure gradients 
in their differential equations are readily found. 
Notice also that normal velocities lie directly on 
the physical boundaries of the solution domain, while 
the tangential velocities and pressure are displaced 
half a cell Interval inside the flovfield. In this 
way the exterior fictitious cells are particularly 
convenient when applying the boundary conditions. 

Figu re 1 also represents the physical problem 
of the first application which contains a downward 
flowing Jet of fluid from a propeller. Initially 
two fluids occupy positions above and below the inter- 
face (called 'therswcllne* ) as shown, so that their 
Initial mass fractions are mj • 1 and •, • 0 (for 
t > Zj, the height of the Interface) and vice verse. 
Also shown Is how the Inlet and outlet flows are 
handled. The exit release flow la via a rectangular 
area In the release structure with the flow rate 





Fig. 1 Fully 3-D Cartesian Crld Schematic. 



Fig. 2 Arrangement of Finite Difference Varlablea In 
a Typical Cell. 

The Bour.lary Conditions 

Finite difference equations simulating the pro- 
blem are set up and solved by way of a time-march pro- 
cess applied to cells within the flow domain of Inter- 
est. Cells touching the boundary thus utilize the 
value on the boundary (In the case of a normal velo- 
city) or values half a cell distant beyond the bound- 
ary (in the case of tangential velocities). 

Interior normal velocity calculations taka the 
zero normal wall values, Che given normal Inlet values, 
on the yet-to-be-dctermlned outlet values as appropri- 
ate boundary conditions during their calculations. 
Interior tangential velocity calculations usa the • 
fictitious values which are placed In the surrounding 
layer of complementary cells. Specification of these 
la after each time-step and after each sweep of the 
cells during the pressure Iteration. With a coarae 
grid, frea-aiip boundary conditions are appropriate 
for tangential velocities, and external values ere set 





tqual to Choir associated Immediately Inter lor valuta. 
On tho other hand, with a (Ina (rid computing through 
tha boundary layer, rto-alip boundary condltlona are 
appropriate for tangential walocltlaa and aatarnal 
valuea are aat equal to tha negative of thalr aeso- 
claced Immediately Interior valuea. 

Specif tcatlona of normal velocltlea at an out - 
flow boundary often poeea a problem, aa It can have 
detrimental upatreaa influence. One might merely 
impose the aero-normal gradient or continual Ive 
condition and aet theae valuea equal to their lmmadl- 
tely upacream valuta (1). When primary Intareat la 
oclng focuaed on the final ateady-atate solution. It 
ftae been found ()) that a auitable constant may be 
added to each auch eatrapulated value, with advantage 
to the rapidity of convergence. Thla constant value 
la chosen ao aa to make the total outlet flus equal 
to the total Inlet flus, thus ensuring the require- 
ment of a macroacroplc mesa balance, (hntlat boundary 
specification la Imposed only after each time-step 
and not after each pass through the ueeh during the 
pressure iteration. 

At planes of symmetry the usual tero normal 
velocity and free-sllp tangential velocity specifi- 
cation are applicable. 

The Inlet value for the turbulent kinetic 
energy k la taken as 0.03u* and the dissipation C la 
calculated from Its definition |e - Cq k ' ' * /t ) and 
the length scale of turbulence t, taken typically 
aa ) per cent 01 the characteristic alse of the 
passage. 

At a plane of symmetry both k and C are given 
the xero-gradlent condition. At a wall, k Is treated 
similarly but c values are specified just Inside the 
computational field at a point P In the terms of the 
local value of k via (A) 
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K.y 
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(5) 


Wall shear stress Is evaluated via the modified log- 
law with the assumption of a constant shear region 
near the wall. Wall functions are employed In u, v, 
w, k and C equations, using the recommended expres- 
sions (A) 


V 

P 


<T/P)„ 





( 6 ) 


E . _JL_ 

(x7pj~ c 172 

w c ll 

where the suffices w and p refer to the conditions 
at the wall and at the nrar-vall grid node respec- 
tively. K and E are constants used in tha log-law 
of the wall. and possess the values 0.A186 and 9.793, 
resepctlvely. 


which la usually of special Interest aa opposed to the 
en route calculations. In Eq. (1) the time-derivatives 
ate approximated by forward one-alded derivatives; 
most apatloal derivatives are approslsMted by central 
differences based on values at t law-level t. Special 
techniques are required In computational fluid dyna- 
mics, however. In the reoreaentatlon of the convection 
terms, and a certain amount of upstream differencing 
la required. 

The difference equations representing Che partial 
differential equations may be written explicitly aai 

♦ - ♦* ♦ At |—1 ( 7 ) 


where f • s, r, w, mj, k and C. Convection, diffusion 
and source terms occur In the parentheses on the right 
hand aide. For the moawntum equations the appropriate 
forms are given In kef. 6. For the turbulence quanti- 
ties k and t, similar terms occur and source term 
Unearitat Ion la effected in the manner recovended for 
always positive variables |«m kef. 7, p. IAS), lm 
all sis of these forward marching aquations, donor 
cell differencing Is used with the convection terms. 

In thla a coefficient a takes a constant value between 
0 and 1, so giving the desired amount of upatreaa dif- 
ferencing (3). A value of 0 gives merely central 
differencing as In tha original MAC code and nuawrlcal 
instability problems arise; a value of 1 glvea tha 
full upstream or donor cell for which, although lass 
accurate. Is atable provided awong ocher thlnga that 
the fluid is not allowed to pass through more than one 
cell in one time-step, see the stability criterion of 
Eq. (11) and other conditions later. 

Although Eq. (S) accomplishes one forward tlma- 
step baaed on conservation of momentum principles, the 
newly calculated velocities will not, In general, 
satisfy the continuity requirement, as expressed by 
the central finite difference form of the continuity 
equation (the first of Eq. (l))l 

S (u ljk‘ t, l-l,j.k > 4 ^ (v ljk' ¥ l,j-l,k > 


4 U (w ljk- ¥ l,j,k-l ) * 0 


( 8 ) 


Terms here are evaluated at time-level t + At. Thla 
Incompressibility condition la imposed by Iteratively 
adjusting the cell pressure. That is, if the diver- 
gence D of a cell la positive (the left hand side of 
Eq. (8) la oosltlve) there la a net mass outflow from 
that cell. This Is corrected by reducing the cell 
pressure. If the divergence Is negative, an Increase 
in cell pressure la appropriate. 

When a cell pressure changes from p to p ♦ Ap, 
the velocity components on the 6 faces of that cell 
change, given from a linear analysis from Eq. (7), 
by the amounts: 

“l.j.k ’ "l.j.k 4 At Ap/(A * P *> 


Solution Procedure 

Finite difference representations are required 
of the governing partial differential equation*. 

The usual Intuitive estimates of one-sided flrat- 
derlvatlves, centered f lrst -der lvat Ives and centered 
second-derivatives are used in representing the 
momentum equations. Superscripts n and (blank) are 
used to denote values at time-level t and t 1 At, 
respectively. Portrayed now are equations enabling 
ona such forward lime-step to be accomplished. Thur, 
atartlng from Initial field values throughout the 
domain of Interest, a time-march process la used so 
aa to advance toward the final steady-state solution. 


“l-l.j.k ’ “l-l.j.k • At Ap/(Ax P «> 

"l.j.k ""l.j.k 4A< Ap/(AjrP,) ( , 

y l,j-l,k ’ V l,j-l,k ‘ At Ap/(A * P «> 

V l.l.k " "l.j.k 4 At * /Wl P «> 

"l.j.k-l * "l.j.k-1 * At * ,(Ab P »> 

Substitution of thesa In Eq. (8) yields tha amount of 
correction to p required as: 
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P - -0,D/|2 At (l/Ax* ♦ 1/Ay* ♦ l/Az*)| (10) 

ul v is the current (nonzero) value of the left 
nan* aide ol [q, (8). Praaaura update Iteration un- 
til tna Da ol all the cells are leaa than cone pre- 
scribed aaall poeltive quantity. The tq. (9) are 
applied with an over-reiaaatlon (actor w between 1 
and 2, a v '.ue near 1.0 often being optlaal, In order 
to speed up the convergence ol the pressure Iteration 
process. After the continuity equation la suffi- 
ciently well satisfied, the values are accepted as 
now-tine values and preparation for the neat atop 
of the tlne-narch can begin. 

Stability and Accuracy 

Accuracy Is enhanced by using snail apatlal and 
tine intervals, at the expense of large conputar 
tine. When total tine requlrenenta necessitate a 
coarse grid. It Is not possible to resolve thin 
boundary layers along confining wall and free-sllp 
BCa tor tangential velocities are sore appropriate 
with linear finite difference appros Inst Ions then 
the expected viscous flow no-slip BCs. 

Once the spatiel subdivision is chosen, the 
tie* lncrenent Bust be restricted in two ways. 
Firstly, it oust be less than (typically equal to 
0.21 to 0.33 tines) the nlnlnun cell transit tine 
taken over all cells: 

At < ■ ln {r* t5t* r} <n> 

Secondly, with a nonzero klnenatlc viscosity, oomen- 
tua oust not diffuse Bore than approximately one cell 
In one tine-step, (or which an estimate (3) Is: 

v At < j ain | Ax 1 , Ay*. As* } (12) 

A similar criterion la given In Ref. 6. 

When the time-step Is so restricted, the re- 
quired amount of ".pstream (donor cell) differencing 
must be achieved oy choosing a slightly larger than 
(typically 1.2 to 1.3 times) the largest of the right 
hand side member of: 

• i • > -« {^. J # i ) <>’> 

where the maximum Is taken over all cells. If a Is 
chosen to be too large, stability Is being achieved 
at the expense of the Introduction of an unnecessari- 
ly larRe amount to dlffusion-llke truncation errors 
(called numerical smoothing^ 

APPLICATIONS AND DISCUSSION 

In order to demonstrate the capability of the 
code Just described, two applications are now given. 
In each case, predictions are given via the use of 
the recommended values (4) of the empirical constants 
In the turbulence model. Typically a rather coaraa 
grid syetea Is used In the demonstrations, although 
clearly finer meshes are needed If accuracy Is 
paramount. Approximately 10 minutes of IBM 370/168 
CP time Is required to allow the final stead, state 
solution to evolve. Results presented are demonstra- 
tlve of the qualitative predictions currently avail- 
able. 

Local Ocstratlf lcatlon of Reservoirs 

The Carton pump consists of a low-energy axial 
flow propeller placed Just below the surface so as to 
provide a downward directed Jet of fluid and thereby 
locally alx reservoirs near the release structure of 


the dam. la thla way high-quality eplllonlon water la 
traneported downwards, so obtaining local deatratlfl- 
catlon and Improved release water quality In the vici- 
nity of low-level release atructures. The flovfleld 
Is fully three-dloens tonal and rhe solution procedure 
neede to include species diffusion and buoyancy forces. 
Figure 3 shows a schematic of the alaing so produced 
when there Is no exit bei earth the propeller. The Jet 



Fig. 3 Scheaatlc of a Typical Propeller Puap and the 
Flovfleld Produced Without Exit Flow. 

or plume penetrates sobs distance below the level of 
the tharaocllne. Figure 4 Illustrates the practical 
application of localized mixing in the proxialty of 
the release structure of a dam with low-level release 
gate. The flovfleld so produced Is fully three-dimen- 
sional . Earlier work by the authors (3) was restrict- 
ed to this axisyaa^tric formulation about a vertlca* 



Fig. 4 Schematic of a Typical Propeller Pump and the 
Flovfleld Produced With Exit Flow Via a Low 
Level Release Structure. 

axis; and good results were portrayed for cases In 
which the propeller was close to the release structure. 
More recently, a three-dimensional version (6) has 
been developed using a constant eddy viscosity aodcl, 
based on round turbulent free Jet theory. Coaparlsoa 
of those predictions with the main source 'Of experi- 
mental data (8) confirmed th.t the main dynamic effects 
are axideled better than the previous two-dimensional 
predictions. The present atudy extends the three- 
dimensional simulation to Include turbulent mixing via 
the two-equation k-c turbulence model. 

For the numerical simulation, a grid ayatem of 
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1*1*1 Internal calls (la a, y and i directions, 
respectively) la employed. (This corresponds to 
lOatall grid ltnee), Ho ajip boundary conditions 
with tha wall (unction approach ara takan at tha 
vertical daai wall and tha horltontal floor o ( tha 
reservoir. fra a slip conditions ara takan at tha 
(raa aurfscs, which la alao assumed to maintain a 
constant height. Tha propallor la posltlonad euffl- 
clantly lar beneath tha surfaca (or aur(aco o((octa 
to bo neglected. Valocllloa aro apocldad a priori 
with (1st prodlas at tha proprllor and salt loca- 
tions. Uaa la aada o( the plana o( aynry through 
tha cantar o( tha propallar and tha oalt| fro# slip 
conditions aro applied thara. Tha available voluoa 
o( tha (low doataln la vary Halted. To allow tha 
outlet o( rolaaar water without draaatlcally decreas- 
ing the (luld level, there la a compensating Inc on - 
Ing (low which la distributed aaong tha top and 
hot too layers In anounts equal to tha eplllonatlc 
ar.4 hypo 1 1 anal lc water released. Tha peripheral 
Ini low is at all elevations In auch a way that 
eplllonatlc water enters Into the upper layer and 
does not disturb tha density profile. The large 
area o( ln(low allows It to be at a low velocity In 
order to avoid any disturbance to the (lowfleld. 
Slallarly, hypollonetlc water la led Into the bo t too 
layer . 

The dilution (actor DF (release water quality) 
la found to ba a (unction of turbulent viscosity |i,, 
denslnetrlc Froude number Frd, nondioensional theneo- 
cline heigh' Zj* , nondlmenslonal (low rata Q*, and 
nondioensional propeller diameter D*, horltontal 
distance K* and depth L*. Results are now discussed 
about the application of tha present 3-0 turbulent 
code to this problem with the following base values 
of the parameters: Frd ■ 2.0, D* - 0.111, K*>0.211, 

L* * 0.211, Q* ■ 2.6, Zj* * 0.6. These values are 
taken so as to represent the same conditions as the 
available experimental data (8). 

The rate of progress of the transient computa- 
tions toward the steady state solution la Illustrated 
In Fig. 3. Prediction of the dilution factor DF as 
a function of Iteration time la shown for both a 
laminar and a turbulent flow with k-C model. The 
latter approach la closer to the experimental evi- 
dence. It may also be noted that steady state con- 
ditions have established themselves In approximately 
17 seconds (corresponding to ii>wi» 320 forward steps). 

In Figure 6 the dilution t ■ >r DF Is shown to 
be a strong function of tha den * trie Froude number 



ri«. 3 Prediction of Dilution Factor DF as a Function 
of a Time ( — — Turbulent, Laminar) 



rig. 6 Dilution Factor DF as a Function of Frouda 

Number Fr^ | Q Experiments ( 8 ), • Predictions 
with k-C Model, 0 Predictions with Constant 
Vlscoalty) . 

Frd In both predlctlva and laboratory data. Tha above 
base values ara ratalned for the other parameters. 

Tha general result la that relaasa water quality 
(fraction of eplllanetic (top) water in tha exit 
stream) Increases with higher values of tha Frouda 
number, such conditions being achieved with hlghar 
Jat velocities from the propeller and/or a lower de- 
gree of stratification in the lake. The figure clearly 
shows the difference between the use of a constant 
eddy viscosity model and the k-C turbulence model! 
the former under-predicting the releesr water quality 
and the latter giving mure accurate predictions when 
compared to the experimental data. From the designer's 
view, the important result is the Froude number at 
vhlch the dilution factor rises sharply, slncie that 
determines the velocity required at the propeller for 
the given conditions. Clearly the 3-D turbulent simu- 
lation is superior to the constant eddy viscosity 
model, and this Froude number is predicted better via 
this approach. The computations were performed with 
the standard k-C turbulence model for the mixing. 

The Influence of buoyancy was via tha buoyancy term 
In the vertical momentum equation. No attempt has 
been made to Include the effects of buoyancy on tha 
generation rates of the turbulence parameters k and C, 
a practice advocated elsewhere (9,10), since there la 
as yet no adequate generalisation of the rather meagre 
evidence. 

Deflected Turbulent Jat 

An Interesting three-dimensional problem arises 
when a turbulent jet enters normally Into a uniform 
steady cross-flow. Figure 7 shows a schematic of tha 
problem, which has previously been discussed (7), 
predicted (11) and experiswntally Investigated (12-13). 
Such problem configurations arise In chimney plumes, 
flow under a V/STOL aircraft, film cooling and dilution 
air mixing in combustor applications. Clearly tha 
problem is mora complicated than corresponding fraa 
jet flows Into quiescent surroundings, which ara axi- 
sytnetrlc parabolic problems. Earlier predictions (11) 
utilised a proprietary three-dimensional Implicit 
computer code to obtain steady state predictions. 

This sophisticated code la more complicated to uaa 
than the present code. Moreover, It la not generally 
available. Hence the capability of the present 
straightforward explicit solution scheme Is now Illu- 
strated for this problem. 






rig. 7 Schematic of the Def lected-Jet Problem (7). 

The simulation la for the turbulent Jet of la- 
let velocity u i , see Fig. 7. Use Is made of the 
vertical symmetry plane. A horizontal top plane Is 
located 18 Jet ladli above the bottom plane. A grid 
system of 7x9x9 Internal cells (la x, y and t direc- 
tions, respectively) Is used. Uniform steady cross- 
flow Is specified at the upstream plane of the stain 
Inflow with velocity u lf> . The u ual tero normal 
gradients are taken at the exit plane, free Blip 
boundary conditions ate used on the top, bottom, and 
side planes. 

Figure 8 shows the predicted velocity vectors 
In the yt-plane through x ■ 0. This gives a clear 
Indication of the behavior, magnitude and direction 
of the deflected Jet. These predictions are for the 
Jet to cross-flow velocity ratio R - uj/u ln • A, and 
tha standard values In the k-c turbulence model are 
used (A). Velocity vector plots of this type allow 
the Jet trajectory (the line Joining points of suixl- 
nua velocity) to be determined. 



experimental work for comparison purposes Inc ludei 
Jordtnson (17), Keffer and Raines (11) |who studied 
the structure of turbulence), Ramsey and Coldeteln (1A) 

| who presented velocity and temperature plots In cross- 
sectional and symmetry planeel, and Chassalng at al (11) 
(who studied both cylindrical and co-axial Jets). Tha 
predicted affect of Jet to cross-flow velocity ratio 
R on Jet trajectories In Illustrated In Figure 9 with 
appropriate experimental data for comparison. Earlier 
predictions (11) are very much In agreement with the 
present study. Considering the scatter In the date, 
the quality of the predictions Is very good, thus con- 
firming that the main dynamic and turbulence mixing 
effects ere modeled adequately in the present cod>.. 



III. 

Fig. 9 Location of the Jet Cental line for Different 
Jet-to-Crossf low Velocity Ratios R “ Uj/u 1# . 

CONCLUSIONS 

A prediction procedure for fully 3-D transient 
turbulent flows has been developed. Based on the Los 
Alamos 2-D SOLA ideas, the transient Reynolds equa- 
tions of an incompressible fluid are solved via their 
associated, finite difference equations directly f.n 
terms of the primitive pressure velocity variables. 
Eddy viscosity Is calculated via the k-C two-equation 
turbulence model. Species diffusion and buoyancy 
forces are Included. The techn 1 ue Is simplified to 
facilitate Its use by persons with little or no ex- 
perience In computational fluid dynamics, and boundary 
conditions, for example, ara particularly simple to 
specify. The code thus represents a basic tool, to 
which user-oriented complexities and sophistications 
can easily be added as required. 


Fig. 8 Velocity Vectors In yt-Plane at x • 0 Showing 
Magnitude and Direction of the Jet (for Velo- 
city Ratio R • Uj/u|„ - 4). 
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Abstract 

P ROBLEMS and progress in the simulation and solution 
of complex turbulent reacting flows for ramjet 
combustor applications are turveyed. The full paper reviews 
the diff. cullies, discusses developments, and demonstrates 
that useful predictions are already being made to aid 
designers. Areas in which further detailed numerical 
prediction research will be most useful include (for the 
nonreacting case) recirculation zone characterization, tur- 
bulence simulation in swirling recirculation flow, irregular 
boundary representation (via analytical or numerical trans- 
formation techniques or finite element methods), and further 
complexities and parameter influences; and (for the reacting 
case) the simulation of complex chemistry, the tur- 
bulence/chemistry/spray interaction problem, and ap- 
plication to realistic three-dimensional problems. 

Contents 

The basic ramiet combustor configuration is the so-called 
sudden expansion dump combustor. In this type of com- 
bustor. liquid fuel is sprayed into the ram air upstream of the 
dump station, although it may also be injected directly into 
the chamber via side wall inlets. Primary flame staoiliza'ion is 
provided by the flow recirculation regions, which may be 
supplemented, at the expense of total pressure loss, with 
mechanical flameholding devices at the air inlct/combdstor 
interface and/or the presence of inlet air swirl, obtained by 
the use of tangential injection ot swirl vanes. The flow 
throughout is multiphase, subsonic, turbulent, and involves 
large-scale corner recirculation zones. With strong rirl in the 
inlet flow a central toroidal recirculation zone (a recirculation 
bubble in the middle of the chamber near the inlet) also 
presents itself. Even gross features of the flow are not known 
quantitatively with certainty: for example, factors affecting 
the existence, size, and shape of the recirculation zones. 1 

The designer has a formidable problem in aerother- 
mochemistry, and the modeling task is to provide a route 
which leads to the accomplishment of design objectives more 
quickly and less expensively than current practice permits. 
Some combustor modeling problems are I) physical 
processes: turbulence, radiation, combustion, and multiphase 
effects; 2) computer programs: 0-, I-, 2-, and 3-D approaches 
ip steady state and transient cases; 3) unresolved problems: 
irregular boundary representation, effect of swirl, recir- 
culation, and wall proximity on turbulence, tur- 
bulence-reaction interaction, multiphase simulation. 

Current combustor design and development problems, the 
needs of the combustion engineer in practice, and proposed 
research tasks which will assist in the attainment of design 
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objectives are becoming clear. Improvements and new 
developments (both experimental and theoretical) can and 
should be made, theoretical modeling being aided by specific 
carefully chosen experiments being performed. 

Mathematical models of steadily increasing realism and 
refinement are being developed, both in the dimensionality of 
the model (together with the computational procedures) and 
in problems associated with the simulation of the physical 
processes occurring. In the modeling and prediction of swirl 
flow combustion, one is involved with simulating the problem 
via simultaneous nonlinear partial differential equations. 
These may be parabolic (boundary layer type) but are more 
often elliptic (tecirculating type) and the solution scheme 
differs according to the category. 

Betv se application of the general partial differ ential 
equali ) is complex, time-consuming, and in a development 
stage, simplified approaches to the problem are extremely 
popular. The most common models include perfectly stirred 
reactors (PSRs), well-stirred reactors (WSRs), and plug flow 
reactors (I'FRs). Model t differ in how these are interrelated to 
Simula'*, -n.ious aspects of the mixing/reaction taking place. 
An important problem in finite rate chemistry is choosing an 
appropriate level of complexity, in view of the large number 
of species and chemical reactions taking place. One solution 
to this problem is the use of the quasiglobal rection scheme 
whose key element is a unidirectional subglobal oxidation 
step. Coupled with this are a number of intermediate 
reversible reactions. The model affords a useful computer 
time saving as compared to a full finite-rate chemical kinetics 
formulation. Well-known models are categorized as integral, 
modular, or hybrid. Modular > - 1 methods give useful 
qualitative trend predictions and, when amalgamated with 
finite difference flowfield predictions via 2-D axisymmetric 4 
or fully 3-D approaches, 5 excellent results are available. 

Most combustion systems exhibit recirculation, and full 
flowfield prediction requires iterative solution trxhnioues. 
Axisymmetric simulations give rise to 2-D elliptic problems 
and involve 2-D storage, stream function-vorticity i-w or 
primitive pressure-velocity p-uv form ilation, and Gauss- 
Scidel point -by-point iteration or line-by-line SIMPLE (semi- 
implicit method for pressure linked equations) iteration 
procedures, the line methods involving the TDMA 
(tridiagonal matrix algorithm). The essential differences 
between the various available computer codes include the 
following: the complexity of the equation set for the 
simulation of the physical processes, the storage 
requirements, the location of variables ir the grid space 
system, the method of deriving the finite-difference equations 
that are incorporated, and the solution technique. In primitive 
pressure-velocity variable formulations a staggered grid 
system is normally used, as recommended by Los Alamos for 
its special attributes. In computational fluid dynamics the 
“best" representation of the convection and diffusion teims 
is essential to the accuracy and convergence of the iteration 
scheme. At high cel! Reynolds numbers a certain degree of 
“upstream differencing” is essential, using, for example, 
these techniques: upwind differencing, a hybrid formulation 
or the Los Alamos zip, donor cell, etc. Sample recent 
numerical prediction studies include Refs. 6 and 7 for the ^-w 
approach and Refs. 8 and 9 for th e p-u-v approach. Excellent 
trade-off between code complexity and quality of flowfield 
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pattern* it ava'lahle via ihe axisyinmetric approach, lor 
example, in Kef. 9 computational tetultt show the interesting 
el fecit of tevetal combustor design parameters, such at 
degree of twirl, performance of a recirculation /one amplifier 
("tup"), and effect of laterally induced secondary air supply 
on subsequent flow field development ar.d combustor pet 
forinance (for example, velocity, temperature, and com- 
position distributions, the occurrence of recirculation /onet, 
and flame ti/e, shape, and combustion intensity). 

Practical ramjet combustor designs exhibit many 
nonaxisymmeiric features mainly because of air and fuel 
inlets at discrete circumferential locations. The requirement lo 
predict a fully I D flowficld arises. Recent significant con 
tnbutioiis to 1 D combustor flow field prediction ate exem 
plified in Kefs. 5 and 10 1 1, where application has been to the 
furnace and gas turbine fields. Computationally, the 
techniques are very similar, except handling very high 
combustion intensities in the latter case causes additional 
problems Implicit tcchnicues for J-D flows with recirculation 
have the advantage of efficiency and stability in their com- 
putations Al! the work cued is implicit in character and tolscs 
directly for pressure and velocities; a staggered mesh system is 
used and some degree of upstream differencing it used for the 
convection terms at high cell Reynolds numbers. Perhaps Ref 
5 is the most ambitious in being a hybrid technique. The finite 
difference solution of some 10 partial uifferential equations 
on a 27x18x7 grid system provides input to a modular 
approach which handles realistic chemical kinetics, with 
presently IJ species undergoing 18 reactions. Encouraging 
results are portrayed for a variety of flow types. 

In deciding and justifying the use of a particular procedure, 
one has some deliberation One of the objectives of the 
piesent full paper was lo clarify the choice and give ap- 
propriate advice, emphasizing computer application where 
appropriate. For example, at the present time the usefulness 
of a fully I D computer code in practice is not clear. Again, 
recent discussions have doubted the industrial value of this, 
especially when current practice and many problem areas 
involve use of 0 1) and I I) models, as successful application 
of modular and hybrid schemes reveal. The 2-D and l-D 
approaches do possess the possibility of being eventually 
capable of a higher degree of realism, but model accuracy is in 
doubt, time and cost requirements are large and, more 
specifically, many current needs do not demand them. 

Model experiments designed to highlight specific sub- 
problems and their inlet actions (and not complex combustor 
Hows) should receive attention lor model development and 
validation. Computational results show some of the in- 
teresting effects of combustor design parameters on sub- 
sequent flowficld development and combustor performance. 
Computational experiments can be, should be, and have been 
performed which complement test cell data describing the 
internal flowficld characteristics Theoretically, it is in the 
area of parameter influences that computer programs show 
their supremacy in terms of time and cost savings as compared 


with experimental work Progress will lead to more realistic 
and cost effective practical combustor modeling This, and 
continued development on the accuracy of (he simulation of 
the physical processes involved, should elevate computer 
modeling to an established place in practical combustor 
design and development progtams. 
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1 . SUMMARY^ OF P ROJ E C T ACT 1 V IT Y_ 

Through the subcontract project number 8178 Dynamics Technology has con- 
tributed to the Oklahoma State University project “Investigation of 
Flowflelds Found In Typical Combustor Geometries" In the following two 
areas: 

(1) Facility assembly In preparation for turbulence measurements 
In a swirling confined jet using a three-wire hot-wire 
probe. 

(21 Co sultution by D. K. McLaughlin on student work at Oklahoma 
r taU 1 1 iverslty. 

2. SUMMARY OF PROGRESS 

Activity 1, 

Considerable effort has been expended in assembling the confined, 
swirling jet facility. Several of the major components of this 
facility: the nozzle, turbulence management section and fan diffuser, 
were constructed by students at Oklahoma State University this past 
summer. Dynamics Technology personnel have been assembling these conv 
ponents on framework constructed for that purpose and mating It to the 
fan supplied by Dynamics Technology. 

A proMem with the diffuser downstream of the fan has demanded our 
attention. The flow exiting the fan has residual swirl. When this 

passes Into the expansion flow In the diffuser a central recirculation 
regions forms. This Is very similar to the flowfleld experienced In the 
test section of the swirling confined jet, which Is the basic phenomenon 
under study. The recirculating flow Is very unsteady and undesirable In 
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the flowffeld upstream of the facility nozzle. 

The recirculation region In the diffuser flow has been eliminated by 
Inserting a section of short straws just downstream of the fan, which 
eliminates the swirl component of the flow. The nonswirling flow then 
separates from the walls of the diffuser as normally expected In 
diffusers with larger than 6* expansion half angle. This separated flow 
has been eliminated by Inserting perforated plates across the flow at 
three locations equally spaced In downstream locations. 

Finishing work on the nozzle and turbulence management sections are pre- 
sently underway. Final facility assembly should take place In early 
March. At that time the swirl pack and probe drive being constructed by 
Oklahoma State University students, will be Incorporated with the main 
facility. 

Planning has begun for the computer software required to control the 
acquisition of data from the three hot-wires on the probe. Analog to 
digital convertors (ADC's) will be controlled with our Tektronlcs 4051 
computer terminal so that numerous batches of digitized signals will be 
Input to computer memory and then transferred to magnetic tape. 
Programs will be written to convert these digitized data Into three 
Instantaneous components of velocity from which statistical estimates of 
mean velocities, and turbulent normal and shear stresses will be 
obtained. Distributions of these quantities will be plotted on the Tek- 
tronlcs 4051 CRT screen from which hard copies will be directly 
obtained. 

The three-wire hot-wire probe requires three channels of anemometer 
electronics which will be supplied by Dynamics Technology. The three- 
wire probe will require extensive calibration at the start, and periodic 
calibration during the series of production experiments. The results of 
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these calibrations are used in the data reduction to convert raw hot- 
wire signals Into the individual velocity components. We anticipate 
beginning the first calibrations In mld-fVjrch. 

Activity 2^ 

Dr. t-tl.aughl 1 n has provid'd consultation services to students at 
Oklahoma State University during the past several months. The bulk of 
this activity has focused on tne MS thesis work of Salim Janjua >nd the 
crossed hot-wire experimental program of Trevor Jackson. This activity 
will continue on an as needed basis. 


